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Abstract. We introduce a class of multidimensional Schrodinger operators 
with elliptic potential which generalize the classical Lame operator to higher 
dimensions. One natural example is the Calogero-Moser operator, others are 
related to the root systems and their deformations. We conjecture that these 
operators are algebraically integrable, which is a proper generalization of the 
finite-gap property of the Lame operator. Using earlier results of Braverman, 
Etingof and Gaitsgory, we prove this under additional assumption of the usual, 
Liouville integrability. In particular, this proves the Chalykh-Veselov conjec- 
ture for the elliptic Calogero— Moser problem for all root systems. We also 
establish algebraic integrability in all known two-dimensional cases. A general 
procedure for calculating the Bloch eigenfunctions is explained. It is worked 
out in detail for two specific examples: one is related to B2 case, another one 
is a certain deformation of the A2 case. In these two cases we also obtain sim- 
ilar results for the discrete versions of these problems, related to the difference 
operators of Macdonald— Ruijsenaars type. 



1. Introduction 

In this paper we consider higher-dimensional analogues of the classical Lame 
operator 

(1.1) L ^ --^ + m{m + l)p{z) , m e Z+ . 

Here p(z) — p(z|1,t) is the Weierstrass p-function with periods l,r. More gener- 
ally, we are interested in multivariable analogues of the so-called elliptic algebro- 
geometric operators L = —(P /dz'^+u{z), which appeared in the finite-gap theory 
initiated in 70's by Novikov j^. This theory provides a beautiful interplay between 
the spectral theory and algebraic geometry, and the Lame operator is the simplest 
and best known member of this family of operators (see for a survey). Since 
then there have been several attempts to generalize some parts of that theory to 
higher dimensions, most notably jSlU, see also |S1E1IZ|. We should stress, however, 
that in general this leads to differential operators with matrix coefficients. On the 
other hand, in |S] it was suggested to consider the quantum elliptic Calogero-Moser 
problem and its versions related to the root systems [5] as natural multidimensional 
analogues of the Lame operator. More specifically, a conjecture from says that 
for integer values of the coupling parameters the corresponding Schrodinger oper- 
ators are algebraically integrable (this is a proper generalization of the properties 
of the algebro-geometric operators to higher dimensions, see ^lE] and Section O 
below) . For the rational and trigonometric versions of the Calogero-Moser problem 
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this was proved in 10 . Elliptic version, however, turned out to be more difficult: 
until now it was known for An case only, due to 

One of the results of the present paper is a proof of that conjecture of ^ for 
all root systems. In fact, our approach applies to a wider class of Schrodinger op- 
erators, which is an elliptic version of the class introduced in see also |13| . 
Their singularities are the second order poles along a set of hyperplanes satisfying 
some special conditions, which encode the triviality of the local monodromy around 
each of the poles. We call the corresponding Schrodinger operators the gener- 
alized Lame operators. The Calogero-Moser operators with integer coupling 
parameters give particular examples of such operators. Our main result says that 
for a generalized Lame operator its algebraic integrability follows from the usual, 
Liouville integrability (in a slightly stronger sense) . The proof uses a criterion from 
based on differential Galois theory. In dimension one we recover in this way 
the main result of m. For the elliptic Calogero-Moser problem the complete in- 
tegrability was proved (for all root systems) by Cherednik 15 , and this allows us 
to prove the conjecture of (Hj. 

A complete description of all generalized Lame operators is an open problem. 
All known (irreducible) examples in dimension > 1 are related to the root systems 
and their deformations which appeared in I16j : we list them all in Section 4. 
We conjecture that they are all algebraically integrable. Using our main result, 
we check this for all two-dimensional examples, since the complete integrability is 
relatively easy to work out in that case. 

One important property of the algebraically integrable operators is that their 
eigenf unctions can be calculated explicitly (at least, in principle). In section 5 
we explain how to find the Bloch eigenfunctions for a given integrable generalized 
Lame operator, in particular, for the elliptic Calogero-Moser problem. As a result, 
we will see that the Bloch solutions are parametrized by the points of an algebraic 
variety, which is a covering of a product of elliptic curves (in a perfect agreement 
with the situation in dimension one, due to Krichever ^;). Let us mention that 
for the elliptic Calogero-Moser problem (in the An case) the Bloch eigenfunctions 
were calculated by Felder and Varchcnko 18 . Our procedure is different and more 
general (at cost of being less effective). We also explain how these Bloch solutions 
can be used to construct the discrete spectrum eigenstates for the Calogero-Moser 
problem. 

In the last three sections of the paper we consider two particular examples of 
the generalized Lame operators in dimension 2, for which we make the formulas for 
the Bloch solutions very explicit. The first example is 

(1.2) i = -A + 2p{x) + 2p(y) + 4p(a; - y) + Ap{x + y) . 

This is a special case of the elliptic Calogero-Moser problem of the i32-type. Our 
second example is 

(1.3) i = -A + 2p{x) + 2(a^ + b'^)p{ax + by) + 2(a + b^)p{ax + by) , 
— l + i\/3sma — — 1 — i\/3sinck - i\/3cosa 



where a is a complex parameter. Such operator was considered by Hietarinta |19| 
who showed that it admits a commuting operator of order 3. Its rational version 
p{z) — z^^ corresponds to a specific choice of the parameters in the family of 
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two-dimensional Schrodinger operators introduced by Berest and Lutsenko j2(J| in 
connection with Huygens' Principle (see section 4 of J3] for details). Note that the 
potential in is real-valued (for real x, y) if a is real and the period r is pure 
imaginary. 

It is more convenient to work with the following 3-dimensional version of (|1.3|) : 

(1.4) L = -df -di-di + 2{al + al)p{aiXi - a2X2) 

+ 2{al + al)p{a2X2 - 03X3) + 2(03 + al)p{a3X3 - aixi) , 

where af + + ^ 0, di ^ d/dxi. Then it is easy to see that L commutes with 
the operator 

Lq = a^^di + 02^82 + a^^dz 
and after restriction to the plane 

a~[^Xi + a2^X2 + a3^"'^a;3 = 

it reduces to the operator with proper a, b. 

We calculate explicitly the Bloch eigenfunctions of the operators 1)1. 2|l . H1.4|l . 
Notice a certain similarity between our approach and the one used by Inozemtsev 
for A2 case |21) . Let us also mention that in dimension two there is a nice theory 
of Schrodinger operators which are finite-gap at a fixed energy level, see |22l I23| . 
It would be interesting to analyze our results from that point of view. 

In the last two sections we also calculate the Bloch solutions for the discrete ver- 
sions of H1.2|) , H1.4|) , which are given by certain difference operators of Macdonald- 
Ruijsenaars type. This raises a natural question about generalizing our results to 
the difference setting. We hope to return to this problem in future. 
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Foundation for Basic Research (grant RFBR-01-01-00803). The work of P.E. was 
partially supported by the NSF grant DMS-9988796, and partly done for the Clay 
Mathematics Institute. 

2. Generalized Lame operators 

Let V — C" be a complex Euclidean space with the scalar product denoted by 
( , ), and A = {a} be a given finite set of afRne-linear functions on V. Let us 
consider a Schrodinger operator 

(2.1) L = -A + u{x), A = d'^/dxl + --- + dydxl, 
with the elliptic potential u of the following form: 

(2.2) w(a:) = ^ c„p(a(rE)|T). 

Here p{z\t) is the Weierstrass p-function with the periods l,r, Im(T) > 0, and 
Cq e C are the parameters which will be specified later (below we will mostly 
suppress r, denoting p(z\t) simply by p{z)). Each of the functions a{x) in standard 
coordinates on V looks as 

a{x) = ao + aixi + • • • + a„a;„ , 
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so a is, effectively, a pair (ao, ao) where aq = Q^(0) G C and ao = grada = 
(ai,...,a„). We assume that each ao is non-isotropic (co)vector, i.e. {ao,ao) = 
flj + • • • + a^^ 7^ 0. Let us project A onto V* , a i-^ ao, denoting by the resuhing 
set. 

We want u to be periodic, more precisely, we assume that the lattice A4 C V* , 
generated over Z by the set has rank < n. For simplicity, let us assume 
that ikAi = n (the general case can be reduced to that by passing to the factor 
space V/AnnM). In that case the lattice C + tC is the period lattice for u, with 
C := Hom(Al,Z) C V. Thus, u may be considered as a meromorphic function on 
a (compact) torus T = V/C + tC, which is isomorphic to the product of n copies 
of the elliptic curve £ = ^1^ + tTj. Singularities of u{x) are the second order poles 
along the following set of the hyperplanes: 

Smg = U U ' ■= • "(^) = ™ + • 

We will assume that all hyperplanes tt™'" are pairwise different; this can always be 
achieved by rearranging the terms in (|2.2(l . This will imply that 

u{x) ^ Ca{a{x) - m - nr)"^ + 0(1) near tt™'" . 

Our next important assumption is that each Cq, in H2.2|l has a form 

Ca — ma{ma + l)(ao, ao) with some TOq, G Z>o . 

Now we are going to put some more restrictions on u demanding its quasi-invariance 
in the following sense. 

Definition. Let us say that the potential H2.2|l as above is quasi-invariant if for 

any hyperplane tt — tt™'" S Sing the (meromorphic) function u{x) — u{sTrx) is 
divisible by {a{x) — m — nr)^™'*"'"^, where s,r denotes the orthogonal reflection with 
respect to n. 

Here are two important examples of such potentials (more examples will appear 
later). 

Example 2.1. Consider the following potential u{x) in C": 

(2.3) u = 2m(m + l)p{xi — xj) , m g Z>o . 

It is invariant under any permutation of the coordinates x^, . . . ,x", and also under 
translations x i-^- x + / for / G Z" + rZ" (Z" is the standard integer lattice in n 
dimensions). As a result, u will be symmetric with respect to any of the hyperplanes 
Xi — Xj = m + nr. Thus, its Laurent expansion in the normal direction will have 
no odd terms at all, hence u is quasi- invariant. The corresponding Schrodinger 
operator H2.1|l is the Hamiltonian of the quantum elliptic Calogero-Moser problem. 
More generally, the elliptic Calogero-Moser problems related to other root systems 
[HI also lead, in the same way, to quasi-invariant potentials. In particular, for the 
rank-one system Ai we have the classical Lame operator 

Example 2.2. This example is in dimension one. The potential u has N poles 
xi,. . . ,xn and looks as 

N 

(2.4) w=:^2p(a;-x0. 
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To ensure its quasi-invariance, one has to impose the condition that u has zero 
derivative at each of its poles, more exphcitly: 

(2.5) "^p'ixi- Xj) ^0 foralH = l,...,iV. 

This system of equations describes the so-called 'elliptic locus' from j^l], which has 
an intimate connection with the classical elliptic Calogero-Moser system and the 
KdV hierarchy. 

Let us call a Schrodinger operator L with quasi-invariant elliptic potential u{x) 
a generalized Lame operator. In trigonometric and rational versions (p(x) = 
sin^^ X or x^^) such operators were considered in |12| . where their eigenfunctions 
were effectively constructed. From the results of the so-called algebraic integra- 
bility of L follows (see the paper for the rational case). We recall the definition 
of the algebraic integrability in the next section, following [H1^2i i^^^ remark 

that in dimension one this coincides with the class of algebro-geometric operators 
which appear in the finite-gap theory, see [3 123 ■ This motivates the following 

Conjecture. The generalized Lame operators are all algebraically integrable. 

As a particular case, this contains a conjecture of (Sj about the algebraic inte- 
grability of the elliptic Calogero-Moser problems. As we already mentioned in the 
introduction, for the A^-case (I2.3|l this has been proved by Braverman, Etingof and 
Gaitsgory in It is also known to be true in dimension one, due to Gesztesy and 
Weikard |14| . In the next section we prove this conjecture under additional assump- 
tion of the usual, Liouville integrability of L (in a slightly stronger sense). As a 
corollary, we will obtain the algebraic integrability of the quantum Calogero-Moser 
problems for integer coupling parameters. 

3. MONODROMY AND ALGEBRAIC INTEGRABILITY 

Let L be a generalized Lame operator as defined previously. Recall that L is 
completely integrable if it is a member of a commutative family of differential 
operators Li — L,L2,...,Ln which are algebraically independent. We assume 
that the Li's have meromorphic coefficients and are periodic with respect to the 
same lattice, which makes them (singular) differential operators on the torus T = 
C"/£ -I- tC The following proposition shows that possible singularities of Li are 
contained in the singular locus Sing of the Schrodinger operator L. 

Proposition 3.1. Let L be a Schrodinger operator regular in an open set U . Then 
any differential operator M on U with meromorphic coefficients commuting with L 
is regular in U . 

To prove the proposition, we will need the following lemma. 

Lemma 3.2. Let S{x,p) be a meromorphic function on T*U — U x V* which is 
polynomial in the momentum p, and {p^,S{x,p)} is a regular function. Then S is 
a regular function. 

Proof. Assume the contrary. The function S* is a finite sum ^ Sk{x)p^, where p^ 
are monomials. Let D C U he the divisor of poles of S. Take a generic point zo of 
this divisor. Near this point D is given by an equation / = 0, where / is analytic 
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at zq and df{zo) ^ 0. Let S = -^(1 + 0(/)) near zq {Q is regular at zq, with 
Q{zo) ^ 0). Then 

{p^ s} = -2fc^K^/-'=-ig + o(r . 

But '^Pi§§-{zo,p) for generic p. Thus, {p^, S} is singular, which is a contra- 
diction. □ 

Now we prove the proposition. Suppose [Af, L] — 0. Assume M is not regular. 
Then we can write M as M' +M" , where M" is regular, and M' has a singular high- 
est symbol. Then [L,M'] = - [L,M"] is regular. Let S{x,p) be the symbol of M' . 
We have {p^, S} is regular (since if it is nonzero, it is the symbol of [L, M']). Then 
the lemma implies that S is regular. This contradiction proves the proposition. 

Notice also that if S{x,p) is the highest symbol of Li, then from [L,Li] = 
it follows that {p^,S} ~ 0. Thus, by Lemma [3.21 S must be regular everywhere. 
Hence, each Li must have constant highest symbol, i.e. Li = Si{d) + . . . for some 
polynomial Si. 

Definition. Let us say that a Schrodinger operator L — —A -|- w in F = C" is 
strongly integrable if the commuting operators Li — L, . . . , Ln have algebraically 
independent homogeneous constant highest symbols si, . . . , s„ and if C[y] is finitely 
generated as a module over the ring generated by si, . . . , s„ or, equivalently, if the 
system si(^) = 0, . . . , s„(^) = has the unique solution ^ = 0. 

Now let L be a strongly integrable generalised Lame operator, so we have the 
operators Li — L, . . . , L^ with meromorphic coefficients on the torus T — C"/£ + 
tC, and Li = Si{d) -t- . . . . First of all, C[V] is locally free as a module over 
C[si, . . . , s„]. This follows from the fact, due to Serre that if f : X Y is a 
finite map of smooth afRne varieties of the same dimension, then 0{X) is a locally 
free O(y)-module. Further, since Si are homogeneous, this module is graded, hence, 
it must be free. Denote by N the rank of this free module. 

Consider now the eigenvalue problem 

(3.1) L^-tJj ^ X.^-tJj , i = l,...,n, 

with A = (Ai, . . . , A„) £ C". Then the space of solutions of this system in any 
simply connected domain in T \ Sing is N-dimensional. So we have a holonomic 
system of rank N onT with singularities along a finite union of hypertori Sing C T. 

Theorem 3.3. The holonomic system (|3.1() has regular singularities. 

Proof. Regularity of singularities is a codimension one condition. Thus, it is suffi- 
cient to restrict our attention to a neighborhood U of a point which lies on exactly 
one subtorus from the pole divisor. Let us assume that this point is 0, and the 
subtorus is locally defined by the equation xi = 0. 

Lemma 3.4. For any r > 0, the singular part of the differential operator {xiYLi 
has the order at most — r — 1, where di is the degree of s^. 

Proof. Let us introduce new variables yj — txj , and write our operators with respect 

to them. Then L = t-^L(°'> + 0(1), t -> 0, where L^°^ = -A+ "'°^'"°+^^ with 

^1 

nia G Z+ corresponding to the chosen hyperplane xi — 0. Let l'"' be the coefficient 
at the lowest power of t in the expression for Li. Then L^"'' is homogeneous, of 
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degree < —di in t (since the symbol of Li has this degree). On the other hand, 
the order of this operator is at most di. Thus, if this degree of L^"' is less than 
di, then its symbol would have to be singular. But Lf^'^ commutes with L'-^^ so 
this contradicts Lemma 13.21 Thus, the degree is exactly di, which implies the 
statement. □ 

Lemma 3.5. Consider a system of differential equations df/dz — A(z)f{z) with 
holomorphic coefhcients on a punctured disk < \z\ < e (/ is a vector function, A 
is a matrix function). Assume that A{z) is meromorphic at z = and that there 
exists an integer-valued function D{i) such that the order of aij(z) at z = is at 
least D{j) — D{i) — 1. Then the system has a regular singularity at z = 0. 

Indeed, let us make a change of variable gi = z^'-'^'^fi. This will change the 
matrix A into a new matrix A which obviously has at most simple pole at z = 0. 

Now let us prove the theorem. First, let us reduce to a first order holo- 

nomic system in a standard way. Choose a collection of N homogeneous poly- 
nomials qi,...,qN in di, which form a basis in C[di, ■■■,d„] as a free module over 
C[si{d), ...,Sn{d)]. Let be a solution the eigenvalue problem (|3.1() . and consider 
the functions 

gi(9)-0,...,gAr(<9)V'- 

Then, for any polynomial q of di, the function qip can be expressed via qiip, with 
the coefficients depending on x and A. Indeed, assume that q is homogeneous of 
degree d. We know that q can be uniquely represented in the form '^gjqj, where 
gj G C[si, s„]. Thus, q'ip — 1j9ji^ij --i Ln)'4' is expressed through differential 
polynomials of ip of degree smaller than d. But gj{Li, ...,Ln)ip = gj{^i, An)?/', 
so, eventually, by induction we get a desired representation of qip. 

Now observe that by Lemma |3.4I the coefficient of qjil> in the expansion of qij; 
has a pole in xi of order at most deg{q) — deg{qj). Thus, we have a holonomic 
system of matrix partial differential equations 

dkqi->P = ^ a[f {x)qjip k^l,...,n, 

which is equivalent to the system (|3.1|l . Each of the matrices a\j'' {k — 1, . . . ,n) 
satisfies the conditions of Lemma [3. 51 with respect to z = xi, for D{i) = — deg{qi). 
In particular, using the lemma, we conclude from the first equation diqiip = ^ ^ij^V' 
that all the solutions have at most power growth in xi when approaching xi = 0, 
hence the system H3.1|l has a regular singularity at xi = 0. □ 

Now take any point xq (z T \ Sing and consider the monodromy of the system 
()3.1(l at point xq, this gives an A^-dimensional representation of the fundamental 
group 7ri(r\ Sing). 

Proposition 3.6. The monodromy group of the system (|3.1|l is commutative for 
any A. 

Proof. Take a hyperplane tt = tt™^" g Sing and let P be a generic point of tt. 
Changing the coordinates if necessary, we may assume that P = (0, . . . , 0) is the 
origin and tt is given by equation xi = 0. From the regularity of singularities it 
follows that there exist 71 , . . . , 7^ G C such that any solution ^ of the system ()3.1|l 
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near P ^ it has a convergent series expansion in the subspace 

r 

0x7^C[[xi,...,a;„]][log(xi)]. 

Substituting such a series into the first equation Lijj = Xiipi one arrives at certain 
recurrence relations from which foUows that (1) we have two 'leading exponents' 
7i = —ma, 72 = ma + 1, both integer; (2) there will be no log(a;i) involved. The 
latter fact is due to the quasi-invariance of u, see |27| for the one-dimensional case 
and J section 2 for a discussion in the multivariable setting. 

As a result, we see that all the solutions are single-valued near tt, so the local 
monodromy corresponding to a loop around tt is trivial. Consequently, the global 
monodromy group will be a homomorphic image of the commutative group ni{T). 

□ 

Corollary 3.7. The differential Galois group of the system (13.1(1 is commutative 
for any A. 

Proof. It is known (see that for a regular holonomic system on a smooth 

projective variety the differential Galois group coincides with the Zariski closure of 
the monodromy group. In our situation the proof is simple: first, we know that all 
solutions of the system are meromorphic in C" . Indeed, this is true outside 
the set of codimension 2, by the proposition above, hence it holds everywhere by 
Hartogs' theorem. Now let F be the field of 'elliptic functions', i.e. meromorphic 
functions on the torus T, and L is the solution field of H3.1|l (on some simply 
connected domain). The monodromy gives the homomorphism 

7ri(r) h-> Dgal{L : F) 

of the fundamental group of the torus T into the differential Galois group. Let G be 
the Zariski closure of the image of this homomorphism inside Dgal{L : F). By the 
main theorem of the differential Galois theory ^29^, to prove that G — Dgal{L : F) 
it suffices to show that in the solution field any G-invariant function is actulaly 
_DgaZ-inavriant, i.e. belongs to F. But any meromorphic G-invariant function is 
7ri(r)-invariant, hence elliptic. Thus, it is Z?ga^-invariant by definition. 

This proves that Dgal{L : F) — G, and the latter is obviously commutative. □ 

Now let us remark (see that a quantum completely integrable system (QCIS) 
on a smooth algebraic variety X of dimension n naturally defines an embedding 

9 : 0(A) ^ D{X), 

where A ~ A" is an affine space and D{X) denotes the ring of differential operators 
on X. More generally, A can be any affine variety with dim(A) = n. Then we 
have an analogous eigenvalue problem 9{g)tp = g{X)tp, \lg G C(A). The dimension 
of the local solution space of this system at generic point of X is called the rank 
of a QCIS. Recall further, that a QCIS S = (A, 6) is algebraically integrable 
if it is dominated by another QCIS S' = {A', 9') of rank one {S is dominated by 
S' if there is a map of algebras h : 0{A) i— s- 0{A') such that 9 = 9' o h). In our 
situation, this implies that apart from the operators Li, . . . , L„, we have additional 
commuting operators which are not algebraic combinations of Li (though, of course, 
they are algebraically dependent with Li). In dimension n — 1 this is equivalent to 
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saying that L is a member of a maximal commutative ring in D{X) of rank one; 
this is known to coincide with the class of algebro-geometric operators. 



Theorem 3.8. Any generalised Lame operator L which is strongly integrable is 
algebraically integrable. 

This follows immediately from the result above and the criterion from . More- 
over, according to (IT] , for generic A the solution space is generated by the quasiperi- 
odic solutions: 

Corollary 3.9. There exist meromorphic 1-forms Uj on the torus T, with first 
order poles and depending analytically on X, such that the functions 'ipj = '^^ give 
a basis of the solution space of (I3.1|l for generic A. 

Each of these functions will be double-Bloch, in terminology of |3(J| : 



(since 4'j has no branching along Sing, these are well defined modulo Z). In Section 
|S1 below we explain how one can calculate these double-Bloch solutions for the 
generalized Lame operators. 

Remark 3.10. Our argument applies to a more general situation when one has 
commuting operators Li, . . . , L„ on an abelian variety, such that the system Liip = 
Xiip, ? = 1, . . . , n, is regular holonomic. Then the triviality of its local monodromy 
around singularities implies that this system is algebraically integrable. 

4. Examples 

4.1. One-dimensional case. Let L be a one-dimensional Schrodinger operator 
L = + u[x), a; e C. Consider the eigenvalue problem Li}) = Xip. Now, 

assuming that u is meromorphic and has a pole at a; = 0, let us consider the local 
monodromy of the solutions of this second order differential equation. Suppose that 
this monodromy is trivial for all A, in other words, all the solutions are meromorphic 
at X = 0. Then, using the classical Frobenius analysis, one can see that u must 
have a pole of the second order, with no residue: u — m{m + -I- 0(1), with 

integer m. Furthermore, in the series for u at x — there must be no terms of 
order 2j — 1 for all j = 1, . . . , to (see 127)). This is exactly the quasi-invariance of 
u with respect to the symmetry x —x. 

Now let u be elliptic, with periods 1, r. Let us demand all the solutions of the 
equation Lij) = Xip to be meromorphic in C (such u are called Picard potentials in 
|25[ll4j '). By the discussion above, this is equivalent to the quasi-invariance of u at 
each pole. More explicitly. 



(3.2) 



ipjix + 1) 



for appropriate aj, bj € V* and for all / G C. Namely, 




M 




mi(mi -|- l)p(x — Xi) , with £ Z. 
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and the poles Xi must satisfy the fohowing system of equations, which generahzes 

(4.1) ^"ijK- + l)p^^''^\x, - Xj) = 

for alH = 1, . . . , M and s — 1, . . . , to^. 

Such L automatically defines a QCIS (of rank 2), and the regularity of singu- 
larities is obvious. Applying the results of the previous section, we conclude that 
L is algebraically integrable. Algebraic integrability of L implies the existence of a 
differential operator P, commuting with L. Since P is not a polynomial of L, we 
may assume that it is of odd order. Thus, L must be algebro-geometric, according 
to the Burchnall-Chaundy-Krichever theory, see e.g. [SI]- Thus, our result in this 
case is equivalent to the main result of ^3] . Notice that our approach easily extends 
to the case of operators of any order (cf. the remark at the end of |14jV 

4.2. Quantum Calogero— Moser system. Let F be a complex Euclidean space, 
dimF — n. Let R — {a} be a reduced irreducible root system in V^*, W be the 
corresponding Weyl group, and the parameters TOq be chosen in a VF-invariant way. 
The corresponding Calogero-Moser operator [H] looks as 

(4.2) L = -A+ ^ Triaima + l){a,a)p{{a,x)) . 

aeR+ 

Let C[V"]^ be the ring of VF-invariant polynomials. By the Chevalley theorem, it is 
freely generated by n elements pi, . . . ,p„, and C[F] is a free module over C[l^]'^, of 
rank \ W\. The following result has been proved in fsl for any PF-invariant iria € C 

Theorem 4.1 (Cherednik). For each homogenious p G C[V]^ there exists a dif- 
ferential operator Lp with the highest symbol p, commuting with L: [L,Lp\ = 0. 
The family {Lp} is commutative. 

For integer TOq, the potential u in H4.2|) is quasi-invariant (cf. Example I2.1|l . 
Altogether this proves the conjecture of |Sj: 

Corollary 4.2. The Calogero-Moser operator (|4.2|l is algebraically integrable for 
integer rua ■ 

Let us mention that for the classical root systems R = A. . . D the complete 
set of commuting operators Li, . . . ,L„ for (|4.2|l was explicitly found by Ochima 
and Sekiguchi [21]. Their results cover the BCn case, too. In that case we have 5 
parameters m, go, gi, g2, gs and the Inozcmtsev operator 

n 3 

(4.3) L = -A-^2m(m+l) ^{p{x^-Xj)+p{xi+Xjj)+^ ^ gs{gs + l)p{xi+uJs) , 

with Ws(.s = 0. . .3) denoting the half-periods 0, l/2,r/2, (1 + t)/2. The Weyl group 
W for this case is generated by the permutations of Xi and sign flips, and Theorem 
14. II still holds true, due to jSl- Applying Theorem 13. 81 we obtain 

Corollary 4.3. The Inozcmtsev operator (|4.3(l is algebraically integrable for any 
integer m, 30,51, 52,53 • 
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4.3. Deformed root systems. Other known examples of the generalized Lame 
operators in dimension > 1 are related to deformed root systems, which appeared in 
Below we describe the set of linear functionals A — {a} and the corresponding 
multiplicities TOq.. 

(1) Ansi'fn) system 13^. 

It consists of the following covectors in C"^"'^: 

Xi — Xj, 1 < i < j < n, with multiplicity (m) , 

Xi — y/mXn+i, i — 1, . . . ,n with multiplicity 1 . 

Here m is an integer parameter, and (m) stands for (m) = max{m, —1 — m}. 

(2) Cn,i{m,l) system 

It consists of the following covectors in C"+^: 

± Xj, 1 < * < j < JT^, with multiplicity {k) , 

2xi, i = 1, . . . , n with multiplicity (m) , 

Xi ± ^/kxn+i, i = 1, . . . , n with multiplicity 1 , 

2^/kxn+i with multiplicity (l) . 

Here k,l,m are integer parameters related as fc = '^21+1 ' ^^^^ i^'^) have the 

same meaning as in A„.i(m) case. In two-dimensional case n = I the first group of 
roots is absent and there is no restriction for k to be integer. 

(3) Here is a _BC„-type generalization of the previous example. Let a-'s(s = . . . 3) 
denote the half-periods, as in BCn case (I4.3|l . The set of linear functionals a E A 
and the corresponding multiplicities look as follows: 

Xi ± Xj, 1 < * < j < with multiplicity (fc) , 

Xi + UJs, i — 1, . . . ,n with multiplicity {nis) , (s = . . . 3) , 

Xi ± ^/kxn+i, i = 1, . . . , n with multiplicity 1 , 

Vkxn+i + ujs with multiplicity (Is) , (s = . . . 3) . 

Here k,ls,ms are nine integer parameters related through k =^ ^2" +1^ * ^ 

... 3. The previous case corresponds to nis = m and Is = I. Again, in case n = 1 
the first group of roots is absent and k may not be integer. 

(4) Hietarinta operator iW\ . 

In this case we have three covectors in C'^, 

a — aixi — a2X2 , f3 = 02X2 — 03X3 , 7 = aixi — 03X3 , 

with rria — mp — m-y — 1. Here are arbitrary complex parameters such that 
al+a^+a^ = 0. Notice that the system is essentially two-dimensional since 
a + P + -f = 0. 

(5) An-i,2{rn) system 

It consists of the following covectors in C"+^ : 

Xi — Xj, ^ < i < j < n, with multiplicity m , 

Xi — ^/rriXn+i, i = 1, . . . ,n with multiplicity 1 , 

Xi — \J —\~ ■mxn+2, « = 1, . . . , n with multiplicity 1 , 

y/rnxn+i — ^/~l — rnxn+2 with multiplicity 1 . 

Notice that for m — 1 this system coincides with the system 2) above. 
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In all these cases a direct check shows that the corresponding potential u is 
quasi-invariant. Notice that in all cases u is symmetric with respect to as soon 
as TOq, > 1. Thus, one have to check the quasi- invariance only for those a where 
iria = 1- 

We believe that all these operators are algebraically integrable. In the next 
section we check this for all two-dimensional examples. 

4.4. Two-dimensional case. Let us check that all known two-dimensional gen- 
eralized Lame operators are algebraically integrable. Apart from the root systems 
^2,-6(72 and G2 considered previously, we have three deformed cases, namely, the 
Ai^i{m) case, the Hietarinta operator and the deformed BC2 case. 
First, let us consider the Ai i{m) case: 

(4.4) L = -dl -dl~dl + 2m{m + l)p(a;i - X2)+ 

2{m + l)p(xi — \/rnx3) + 2{m + l)p(a;2 — V^^s) ■ 

In this case we can use the result of 'S^' where the complete integrability of (14.41) 
was established. First, L obviously commutes with Lq = di + 82 + -^^s. 

Proposition 4.4 For any m there exists a third order operator L2 commut- 

ing with Lq, L, and its highest symbol is + 82 + y/md^. 

It is easy to check that as soon as m 7^ — 1, the highest symbols of Lq, L, L2 will 
satisfy the requirements of the strong integrabilty (notice that for m = — 1 the op- 
erator L is trivial). As a result, for integer m we obtain the algebraic integrability 
of the operator (|4.4|l . Note that for the special case m — 2 the algebraic integra- 
bility of H4.4|l was demonstrated in j33J by presenting an explicit extra operator 
commuting with Lq, L, L2. 

Now let us consider the Hietarinta operator 

(4.5) L = -dl - dl -dl + 2[a\ + a\)p{a^x^ - 022:2) + 

2(02 -I- a\)p{a2X2 - 03X3) -f 2(03 -I- a\)p{azxz - oiXi) , 

where a\ + ak, = 0. Such L commutes with Lq — a'{^d\ -|- a2^d2 + ag ^^3, and 
we need one more operator for the complete integrability (we assume that all ai are 
nonzero, otherwise L is reducible). Such operator was found in |19| . 

Proposition 4.5 (|19p. There exists a third order operator L2 commuting with 
Lq, L above, and its highest symbol is aidf + 0282 + 0,38^. 

If we denote the highest symbols of these three operators as si, S2, sa, then it is 
easy to check that the system si(^) = S2{£,) — S3(^) = has the only solution ^ = 0; 
the only exception is the case when 0^ = 02 = 03. As a result, we conclude that 
for all values of the parameters (apart from the case of = = flg) the Hietarinta 
operator is strongly integrable, thus, it is algebraically integrable. Note that this 
also follows from IS' where one more operator commuting with Lq, L, L2 was found. 

Remark 4.6. In the case a\ = = a\ the operator (|4.5() is still algebraically 
integrable, although it is no longer strongly integrable. 

Finally, let us consider the deformed BC2 case. The Schrodinger operator L has 
the following form: 

(4.6) L^-8l-8l + U{x,y), 
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where U = 2{k + l){p{x + \/ky) + p{x — \/ky)) + v{x) + w{y) and v, w are given by 
the expressions 

(4.7) w = ^TOs(™s + l)p(a; + tJs) , w = ky^^lsjls + l)p{-\/kv + ujs) ■ 

s s 

Here wq, a;i,a;2, ^3 are the half-periods and ms,ls and k are nine parameters such 
that k = (2ms + l)/(2/s + 1) for all s = 0, 1, 2,3 (thus, effectively, L contains five 
independent parameters) . 

Proposition 4.7. For any values of the parameters rUg, Is, k such that k = ^2" 
the following operator commutes with L: 

M^^dl- kd^ + 2{U - w)dl - iVkik + l){p{x + Vky) - p{x - Vky))d,dy+ 

2k{U - v)d^j + {2v' + 2{k + 1)(2 - k){p'{x + y/ky) + p'{x - Vky)))d^ + 

{2k^w' + 2^/k{k + l)(2fc - l)(p'(x - Vfcy) - p'{x + Vky)))dy- 

(fc + i fpix + Vky)pix - Vky) + 8(P + l)(p^(a; + Vky) + p^{x - Vky))- 

A{k + l){v + kw){p{x + Vky) + p{x - Vky)) + v" - + k{w" - w^) . 

Here v',w' and so on are the derivatives with respect to the corresponding variable. 

Now since the system ■Ci + = + ^^1 = ^ does not have nontrivial solutions 
as soon as k ^ —1, we conclude that L is strongly integrable (for = — 1 this is 
also true because L is reducible in that case). Hence, for any integer values of the 
parameters ls,ms the operator l|4.6|l - 14.7|l is algebraically integrable. 

5. Bloch solutions 

Let L be a generalized Lame operator which is strongly integrable, thus alge- 
braically integrable. We know already that for generic A the solution space of 
(j3.1(l is spanned by the meromorphic double-Bloch solutions. Now we are going to 
explain how one can, in principle, calculate them. 

Let yV denote the following linear subspace in the space of meromorphic functions 
on C". First, its elements are holoniorphic everywhere apart from the singular 
locus Sing — [J tt™'" of the operator L, where they may have poles, of order < 
along TT™'". Next, take any hyperplane tt = tt™'" with Stt denoting the orthogonal 
reflection with respect to tt. Then any function ip Cz W must have the following 
property: 

(5.1) 93(0;) - (-l)™°^(s^a;) is divisible by (a(a;) - m - nT)™°+^ 

Proposition 5.1. The subspace W is stable under the action of L: i(yV) C W. 
Furthermore, any meromorphic eigenfunction ip of L must belong to W. The same 
is true for any of the commuting operators L2, . . . , Ln: Li{W) C W. 

Proof. Consider any hyperplane tt G Sing and adjust the coordinates in such a way 
that TT is given by equation x„ = 0. Take a generic point in tt and expand ip in Lau- 
rent series in normal direction to tt, i.e. ip = J^jez'^ji^nY , a,j = aj(xi, . . . ,a;„_i). 
Now put 

M := {-TO„ -I- 2Z>o} U {m„ + 1 + 2Z>o} 

and split the series into two parts, with j € M and with j e Z\M: ip = ipi + ip2. 
First claim is that an application of L to (pi will produce a series of a similar kind. 
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This follows directly from the quasi-invariance of u, proving the first part of the 
proposition. On the other hand, if aj(a;„)-' is the first nonzero term in if 2, then 
applying L to ■!/'2 will give a series starting from (a;„)-'^^, which would contradict 
the equation Lip = Xip^ thus proving the second claim. 

In a similar way, if L'L = LL' for some other operator L' , then L'{LY = {LY'L' 
for any r > 1. Thus, if W := L'{W) then L^'iW) C W' for any r. Now suppose 
we could find a function ip in W" which is not in W. Then, consider a series 
expansion of p in the direction, normal to tt, and split it into </? = (/3i + (/32 as above. 
If ip2 7^ 0, then L^ip would have a pole of an arbitrarily high order along n (as 
r increases), which is impossible for an element in W" (since W" = L'W and L' 
has meromorphic coefficients). This would contradict the inclusion L^{W') C W' . 
Thus, W = W. □ 

Now let -0 be a double-Bloch solution of H3.1|l . so for appropriate a,b G V* and 
for any I £ C we have: 

(5.2) ib{x + = V(a;)e2"*<°''> , ^{x + tI) ^ ^{x)e^^''^^^^^ . 

We know that tp is meromorphic in C" with possible poles along the hyperplanes 
TT™'", of order m^. In our discussion below we restrict ourselves to the case when all 
the linear functions a & A have zero constant term, i.e. a = £ V* for all a, so 
a{x) = {a,x). Everything extends to the general case with obvious modifications. 
As a result, we see that ip can be presented in the form 

(5.3) 5{x)=Y[e{{a,x)r'-, 

for some holomorphic <I>(x). Here 6* = 6*1 is the classical (odd) Jacobi theta function. 



(5.4) e{z) = Y^ exp(7ri(n + 1 /2) V + 2Tri{n + 1/2) (0 + 1 /2)) . 

Recall that 6{z) has the following translation properties in z: 

e{z + 1) = -e{z) , 9{z + t) = ~e-'^''"''''''9{z) . 

This determines the translation properties of S. To write them down, it is convenient 
to introduce the following linear map Q : V V* which is defined as 

(5.5) fl:x^Y^ ma{a, x)a . 

Note that Vl maps the lattice £ to a sublattice in = IIom(£,Z). We also need 
a covector Q = \ '^aeA "^qQ^- 

Under these notations, we have the following translation formulas for any I e C: 

(5.6) ,5(a; + /) = e2"<e,/>j(2.)^ 

(5.7) 6{x + It) = e2'^''(e,0-2^''(™,^>-^»(™,0^j(2.) 

As a corollary of (|5.2(l and (|5.6(I - H5.7|I . we conclude that the numerator $ in (|5.3|l 
must have the translation properties as follows: 

$(x + /) = e2"<'^+^?'')$(a;). 
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The vector space of entire functions with such properties is finite-dimensionaL In- 
deed, let w{x, y) denote the bilinear form on V associated with the operator 17: 



(5.8) 



^{x,y) ^ (a, x) {a,y) . 



It is symmetric positive and integer- valued on the lattice C. Let us define B[^] by 
the following series: 



(5.9) 



6 



(x) 



lec 



exp(27ria;(Z + p, x + q) + TTiTuj{l + p,l + p)) . 



It has the following translation properties: 



(5.10) 



(5.11) 



e 



(x + o = e2""('^p)e 



(^) 



e 



{x + It) — e-2iriw(;,a:+ij)-7rirw(;,i)Q 



(x). 



It is easy to show (see e.g.|3S]) that the space of holomorphic functions with such 
translation properties has dimension equal to [A^ : ^C]; this is equal to det(riij) 
where fiy = {ilei,ej) for some basis ei, . . . , e„ of C. A natural basis in this space 
is given by the functions 6[q^J with r e il^^Ai running over the set of represen- 
tatives in U~^{M)/C. 

For the later purposes, let us use slightly different basis, namely, the functions 



(5.12) 



e<'^''^>e(x + 7 + r) 



e := e 



ren-\M)/C. 



It is easy to relate the parameters k ^ V* and j ^ V to p, q: 

k — 2'niQ.p , 7 = q + . 

Let us denote the linear space generated by the functions H5.12|l as Uk,-!- 

We conclude that -0 must belong to linear space 6~^Uk,'f with fc,7 related in 
a simple way to the 'quasimomenta' a, 6 in H5.2|) . Now recall Proposition 15.11 It 
implies that ip must belong to the (finite-dimensional) subspace Wk.-y = [6^^Uk,'y) H 
W. It also implies that L (as well as any of i^'s) preserves this finite-dimensional 
space, so we can eventually find ip by diagonalizing the action of L on Wk.-y Note 
that since the double-Bloch solutions ■0 form an n-parametric family (with A = 
(Ai, . . . , A„) being the parameters), this space Wk,^ will be nonzero only for (fc, 7) 
belonging to a certain rt-dimensional subvariety. In most cases dimWi;.^ < 1, so 
the (unique) function ipk.-y generating Wk^-y 7^ will be an eigenfunction for L 
automatically. 

All this simplifies a little for the Calogero-Moser models, so let us consider this 
case in more detail. For a given reduced, irreducible root system R in V* and a 
fixed VF- invariant m^a), we consider the Calogero-Moser operator 



(5.13) 



L : 



-A- 



rriaima + 1)(q!, a)p((a, x)|t) 



The Bloch solutions must be of the form ip = S ^<I> as in H5.3|l . Note that 6 in this 
case has the following symmetry: 

(5.14) S{wx) = em{w)S{x) foranyweVF, 
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where Sm is the one-dimensional character of W such that 

(5.15) e,n{sc.) - (-1)"° . 

The biUnear symmetric form (|5.8|) is obviously VF-invariant, and since R is ir- 
reducible, uj must be proportional to the W-invariant scalar product {x,y). So, 
Lu{x,y) = K ■ {x,y) for some k which depends on rUa. (For instance, if R consists 
of one W-orbit only and rria = m, one has u}{x^y) = mh{x\y) where h — h(R) 
is the Coxeter number and the form {x\y) on V is normalized in such a way that 
(a^|a^) = 2for all S i?^.) 

The lattice £ in this case is the coweight lattice of R, while A4 is the root 
lattice Q. Still, the numerator $ must belong to the finite-dimensional space Uk,-y 
spanned by the functions H5.12|l . 

We already know that a Bloch solution tp = (5^^$ appears only for those k,j 
when S~^Uk,f n W 7^ 0. Such A:, 7 can be effectively determined. Indeed, due to 
(|5.14(l . the conditions H5.1|l for tp reduce to the quasi-invariance of <&: 

(5.16) ^x) - <^{sax) is divisible by (a, xf"'"+'^ for aU a e i? . 

(These are local conditions near tTq, — {x : {a,x) — 0}, similar conditions for other 
hyperplanes tt™'" G Sing will follow because tp is quasiperiodic.) 
These conditions can be rewritten as 

(5.17) (a^, 9)2^-1$ = for(a,2:)=0 and j = 1, . . . , m„ , 

with {a^,d) denoting the derivative in a^-direction. 

Now recall that we have the period lattice C = , and $ belongs to the linear 
space Uk,f of the functions with the translation properties (|5.10(I ~ H5.11() . Consider 
the following sublattice C" C C: 

(5.18) /:":=/:'e/:", C' = Cnkeia, C" = CnRa'^. 

Let hl^^ denote the space of theta functions with the same translation properties, 
but for the translations I from C" only. Obviously, we have a natural inclusion map 
^fe,7 ^ ^fc 7- It is possible to describe this linear map explicitly, using the standard 
bases in both spaces (look at the formula (|6.4|l below which is a particular example 
of such relation). According to (|5.18(l . the lattice is the direct orthogonal sum 
of two sublattices. Thus, the corresponding theta functions from U^^ will be the 
products of the {n — l)-dimensional theta functions related to £' and the one- 
dimensional theta functions related to the lattice C" . This corresponds to the 
decomposition oiU^^ into a tensor product: U^^ — iU^^)' ® iU^^)" . Applying 
derivative in oi^ direction will affect the one-dimensional theta functions only. As 
a result, for each j we have an expHcit linear map F"'-' from Uk 7 to {U^^)' given 

by 

F"'J"^=(«\9)2J-Vl(„,.>=o- 
This map is given by a matrix whose entries are certain combinations of one- 
dimensional theta- functions and their derivatives. Now we can organize all these 
maps for a G P,+ into one big linear map 

(5.19) r-.Uk,^^ (U^^^Y, r = (F"-^)„efl^,,=i,...,™„, 

j = l...ma, 

with F"'-' defined above. We can think of F as an M x matrix, where A^, M are 
the dimensions of the source and the target spaces, respectively. 
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The outcome is the fohowing: a double-Bloch solution -0 appears exactly for 
those fc, 7 where this linear map has nontrivial kernel. This gives equations on such 
A:, 7 (by equating to zero all x iV minors of the matrix T). In its turn, the kernel 
will determine a corresponding Bloch eigenfunction. (If the kernel has dimension 
> 1, it still defines an invariant subspace for the action of i, so we have at least 
one double-Bloch solution.) 

So, let C denote an analytic subvariety in V"* x given by 



Formulas H5.10|l - (|5.12|l make clear that C is invariant under the following transfor- 
mations of fc, 7: 



Now, for a hmction of one variable f{z) — e''^9{z + c) its derivatives at 2; = are 
obviously polynomial in k. Therefore C, after being factored by the translations 
above, can be considered as an algebraic covering of an abelian variety (a product 
of elliptic curves) + t£' where £' = Q-^M. 

Proposition 5.2. The double-Bloch eigenfunctions of the Calogero-Moser operator 
(|5.13|l are parametrized by the points of an algebraic variety which is a covering of 
an abelian variety €.""/£' + tC' where C! = fJ^^Q, Q is the root lattice and the map 
f2 is defined by (|5.5|) . 

A similar analysis applies to any (integrable) generalized Lame operator, so the 
double-Bloch eigenfunctions are also parametrized by the points of an algebraic 
variety covering a product of elliptic curves. 

Now let C denote the result of factoring the variety H5.2()|l by the translations 
(j5.21|l - H5.22|l . It is an algebraic variety parametrizing the double-Bloch eigenfunc- 
tions of L. Below, following JSj, we will refer to it as the Hermite— Bloch variety 
for L. It differs from the complex Bloch variety, traditionally defined as the set of 
(/I, £■) e (C^ )" X C such that there exists ^ with Lt/j = E^p and tp{x + li) ^ t'-ii'i^) 
where /i, . . .Z„ is a basis in C. Note that the latter is a transcendental complex 
analytic variety. 

Remark 5.3. In case of the Calogero-Moser operator related to a root system R, 
there is a natural action of the Weyl group W on the Hermite-Bloch variety C. 
Also, there is a natural projection of C onto C" sending = Vfc,7 to the set of 
eigenvalues Ai, Litp = Xiip. This is a |W^|-sheeted covering and the Weyl group acts 
on C by permuting the points in the fiber. 

Remark 5.4. Note that our results do not contradict the theorem of Feldman- 
Knorrer-Trubowitz j36| in dimension two, since their result only applies to a real- 
valued smooth potential m in M^. 

The Hermite-Bloch variety C is a subvariety in the total space of a certain bundle 
over the product of elliptic curves defined by H5.21|l - H5.22|l . This bundle naturally 
compactifies to a bundle with the fibers isomorphic to the projective space P". As a 
result, C compactifies to a projective variety, covering the product of elliptic curves. 



(5.20) 



C = {(fc,7)| kerr^O}. 



(5.21) 
(5.22) 



(fc,7) ^-^ (fc,7 + /), len-^M, 
{k,j)^ {k + 2TTini,-f + Tl), lefl^^M. 
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The variety C is closely related to the so-called spectral variety which is defined 
as follows. Suppose L is a strongly integrable generalized Lame operator, so we have 
n commuting operators Li = L, . . . , Ln, which generate a commutative subalgebra 
in the ring of PDO with meromorphic elliptic coefficients. Then by 11 , theorem 
2.2, the centralizer of this subalgebra will be a maximal commutative ring which 
we will denote by Z{L) (using Proposition 15.11 one can show that the operators in 
this ring will share a common family of the double-Bloch eigenfunctions). Each 
operator in Z(L) must have constant highest symbols, by Lemma 13.21 Then from 
the strong integrability we immediately derive that Z{L) is finitely generated. Thus, 
SpecZ(L) defines an affine algebraic variety, which we call the spectral variety. It 
is not quite clear whether the spectral variety is isomorphic to the Hermite-Bloch 
variety (for instance, the latter may not be affine), but at least they must be 
birationally equivalent. 

Finally, let us remark on some algebraic geometry behind the double-Bloch so- 
lutions and the Hermite-Bloch variety for the Calogero-Moser system (|5.13|) . We 
consider the torus T = C"/Q^ + 'tQ'^ where = i?^ ® Z is the coroot lattice. 
Let us define the following subsheaf Q C 0{T) of the structure sheaf of T by re- 
quiring its local sections to have zero normal derivatives of order 1,3,..., 2toq — 1 
along each of the hyperplanes {a,x) =0 (considered as hypertori in T). The sheaf 
Q can be considered as the structure sheaf 0{X) of a singular variety X, with T 
being its injective normalization (cf. [SH). Such X is projective; it is a |VF|-sheeted 
covering of the weighted projective space T/W considered by Looijenga [37], see 
also Notice that from the results of 00] it follows that X is Cohen-Macaulay 
and Gorenstein. Let us consider now the group Pic(X) of invertible sheaves on X. 
Then each of the double-Bloch solutions ipk,^ represents a meromorphic section of a 
degree zero line bundle on X (to define degree, we use the pull-back to the torus T). 
In this way the Hermite-Bloch variety for the Calogero-Moser system becomes an 
n-dimensional subvariety in Pic^[X). It would be interesting to study this relation 
in more detail. 

Remark 5.5. An interesting thing is to analyse how the spectral variety changes 
when T goes to +ioo (trigonometric limit). In this limit the spectral variety becomes 
rational and is relatively well understood. Thus, one could think of the whole family 
depending on r as a deformation of this rational variety. This point of view was 
used in (41) to construct the spectral surface in the simplest A2 case. 

5.1. Discrete spectrum eigenstates. Let us explain how the Bloch solutions 
can be used to construct the discrete spectrum eigenstates of L. Our discussion is 
strictly confined to the Calogero-Moser operator (|5.13() . We take a purely imaginary 
T, this ensures that the potential in H5.13|l is real- valued for x G M". The Calogero- 
Moser operator L is defined on a dense subset of i^(M") and it is self-adjoint only 
formally, and its Bloch solutions are singular. It has square-integrable eigenstates, 
though. Namely, let ^ = ^k,-y be one of the double-Bloch solutions constructed in 
the previous section. Given such a ip, let us symmetrize it as follows: 



where W is the Weyl group of the root system R and Em is the character (|5.15|l . The 
Calogero-Moser operator L is W^-invariant, thus the constructed ^E* will be again 
its eigenfunction (by the same reason, it will be an eigenfunction for all commuting 



(5.23) 
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operators Li). A priori, ^E* might have poles in R" along the hyperplanes (a, x) = c, 
c E Z. However, it is easy to see that 4' has no poles along the hyperplanes 
(a, a;) = 0. This follows immediately from the properties (|5.1|l of ip. To avoid the 
appearance of singularities on other hyperplanes (a, x) = c, one has to impose the 
condition that all the terms ip{wx) in the sum H5.23|l have the same Bloch-Floquet 
multipliers with respect to a shift x x + I with I € C = P^. This means that 
exp{k,l) = exp(wfc,Z) for all w € W and / £ , which in its turn implies that k 
belongs to the lattice 2TiiP. So, we have the following result. 

Proposition 5.6. Let L be the Calogero-Moser operator H5.13|l . Then for any 
point (fc,7) of its Bloch-Hermite variety which satisfies an additional condition k G 
2'KiP ( with P being the weight lattice for R ), the corresponding function H5.23|l ( if 
nonzero) will be a nonsingular in M" eigenfunction of the Calogero-Moser operator 
(|5.13|l and of the higher operators L2, . . . , Ln- 

By construction, \1/ vanishes along the hyperplanes (a, x) G Z, and it gets a factor 
of (— under the orthogonal reflection with respect to such a hyperplane. 
Since these are the reflection hyperplanes of the affine Weyl group of R, they cut 
R" into its fundamental domains (alcoves), so the restriction of \1/ to each alcove 
will be, essentially, the same. We can restrict ^ to any alcove, extending it by zero 
outside, and this gives us a finitely supported smooth eigenfunction of L (notice 
that in the complex domain it still has poles). We see from this that the discrete 
spectrum of L in L^(R") is infinitely degenerate (one says that the spectral problem 
for L splits into identical spectral problems on each of the alcoves) . Morally, this is 
the reason why one should expect the same spectrum considering L not on L^(R") 
but on the space L^(T)^ of PF-invariant functions on the torus T = R"/(5^, as in 
. The latter case is simpler from the technical point of view, since the operator 
L is essentially self-adjoint on L'^(T)^ , see 0^1 for the details. 

In Komori and Takemura considered the elliptic Calogero-Moser problems 
as a perturbation (in t) of the trigonometric case r = +i(X), and Theorem 3.7 of 
[IB] claims that for sufficiently small p = e^^^'^ the family of eigenfunctions (Jack 
polynomials) which corresponds to p = 0, admits analytic continuation in p, and the 
resulting functions will give rise to a complete orthogonal family of eigenfunctions 
of L in L^{T)^ . One can show that our family in the limit r +ioo specializes to 
the Jack polynomials. Comparing this with the previous discussion, we conclude 
that our family must coincide with the one considered in |46|. 

6. Calogero-Moser model of B2 type 
In this section we consider the following 2-dimensional Schrodinger operator 

(6.1) L = -A + 2p{xi) + 2p{x2) + 4p(xi - X2) + 4p(xi + 2:2) , 

where p(z) = p{z\t) is the Weierstrass p-function with the periods 1, t (Imr > 0). 
Our goal is to calculate its double-Bloch eigenfunctions, i.e. such ip that 

(6.2) -ipix + Cj) = Xjijj{x) 

(6.3) V(a; + Tej) =A^,^(x) (i-1,2), 

where (ei, 62) is the standard basis in and (Ai, A2, /ii, /i2) are fixed Bloch-Floquet 
multipliers. 
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First we recall some standard definitions and formulas from the theory of theta- 
functions, see |43l 135) . Let 6'[^] be the one-dimensional theta-function (with char- 
acteristics), defined by the following series: 



e 



(z|t) — ^ exp(7ri(7i -f a)'^T -f 2iii{n + a){z + (3)) . 



Notice that a and (3 are defined modulo 1: 



"a + I 


= 9 


a 




a 




a 


P 




P. 




_/?+ 1_ 




P. 



Later we will need the following formula which can be easily derived from the 
definitions: 



(6.4) 



{xi\t)9 



a2 




(^2|r) 



a+ 




{x+\2T)e 

Ck+ + 








{x^\2t) 
{x+\2t] 



(x_|2r), 



where a± — \ {ai ±0:2), x± — xi ± X2- 

We will mostly use dlY^'^ which we will denote simply by 0{z), which will always 
stand for the odd Jacobi theta function H5.4(l . 

According to the previous section, ijj must have the form 

<^{X1,X2) 



(6.5) 



eixi\T)e{x2\T)9{xi - X2\T)e{xi + X2\t) 

Here $ is nonsingular in C^. The translation properties for tp easily translate into 
the properties of <J>: 

$(x + e,) = -Xj^x) , 

^{x + TCj) = -/ije-3'^'^-6"^^$(x) . 

Standard considerations from the theory of thcta-functions show that the linear 
space of functions with these properties has dimension 9 and $ must be of the form 



i/i 




(3x1 + 7i|3t) 



j73- 





(3.T2 + 72|3t) , 



(6.6) = el^v{KlXl + K2X2) ^ c,y6» 

0<i j"<2 

where Cij are arbitrary constants and parameters 7^ , Kj relate to Xj , Hj as follows: 



(6.7) 



A, 



Remark 6.1. The shifting 

(6.8) 7, ^ 7, + 1 

does not change the space (|6.6f) . the same applies to the shifts 

(6.9) {K,,j,)^{K,+2m,j,+T). 

Conversely, for any given Xj , /ij the corresponding (7^ , Kj ) are determined uniquely 
modulo shifts 
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Now, in accordance with Proposition 15.11 we impose the foUowing 'vanishing' 
conditions on <i>: 



(6.10) 
(6.11) 
(6.12) 
(6.13) 



ai$ = forxi^O, 

92$ = for ^2 = , 

(^1 + ^2)$ = for + a;2 = , 

{di - ^2)$ = for xi - = . 



As we will see below, for a certain 2-diniensional surface in 4-diniensional space of 
parameters {kj,aj), the conditions (|6.1U|I - (|6.13|I cut a one-dimensional subspace in 
9-dimensional space (|6.6|l . Thus, the corresponding -0 will be an eigenfunction for 
L automatically. 

To determine the corresponding [Kj , 7^ ) , let us rewrite $ using (|6.4|l and making 
identification 6'["+i] = 0['^]: 



0<i/m<2 

7/3+1/2' 






7/3" 








(3a;^ 



(3a:+ +7+|6r 
f 7+|6r 



m/i 




to/3 + 1/2' 




(3a;_ + 7_ |6r) 
(32:_ + 7_ |6r) 



where K± = ^{Ki ± 7± = 71 ± 72 and 

(6.14) c/„j = Cij with i = I + m (mod 3) , j = I — m (mod 3) . 

It is easy to see now that (|6.12|l leads to six linear equations on c/„i of the form 



1=0 



iCln 



0, 



E^' 







(m = 0,l,2) 



for certain explicitly given A = {Ai), B = (Bi). For generic parameters 71,72 the 
vectors A, B will be linearly independent. Therefore, these 6 equations determine 
the 2-dimensional kernel of the matrix C — (q„i). Similarly, H6.13|l gives six more 
equations for q™, which determine the cokernel of C. Thus, for any Ki,K2 and 
generic 71,72 the vanishing conditions Hfi.l2|l - H6.13(l determine cim (and, hence, $) 
uniquely up to a common factor. In principle, it is straightforward to write down 
explicit expressions for the coefficients but they are cumbersome and not very 
useful. However, there is a better way of getting an expression for by taking a 
limit in the formula for the difference case, see (|7.7() below. It turns out that 
<& has the form 



(6.15) 



eia,\r)eia2\r) J- 



where the parameters fci, k2 and oi, 02 relate to Kj, jj as 



(6.16) 



7,-(l + r)/2. 



The coefficients 6^ depend on x and ai, a2 and are given by the following recipe. Let 
us introduce formal commutative variables A, B, C, D. We also need the following 
scalar A depending on r: 

A :=0"'(O|r)/0'(O|r). 
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Now introduce U,V a,s U -.^ A + B - C , V -.^ A + B - D and put 

522-1, &21=2[/, bi2^2V, 

(6.17) b2o = -X + U\ bo2 = -X + V\, bn = WV , 
bw=2V{-X + U^), boi = 2U{-X + V^), 
6oo-(-A + C/2)(-A + y2). 

After that one opens the brackets, so each of bij becomes a sum of monomials 
in A, B, C, D with scalar coefficients, and then replaces each monomial using the 
foUowing rule: 

(6.18) cAPB^C'D' — > c6l(P) (xi + ai)6'(«) {x2 + 02)6'^''^ (xi - X2)e^'^ (xi + X2) . 

Here 9 — 6'(z|r), as before, denotes the odd Jacobi theta function, and the upper 
index in brackets refers to taking derivatives in z. We treat a scalar as a multiple of 
A'^B'^C'^D^ assuming, as usual, that /'^''^ = /. To iUustrate this, we present below 
some first of the coefficients: 

622 =0{xi + ai)9{x2 + a2)6{xi - X2)9{xi + X2) , 

^21 =26' [xi + ai)e{x2 + a2)e{xi - X2)B(xi + 2:2)- 

(6.19) 2e{xi + ai)e'{x2 + a2)e{xi - X2)e{xi + X2)- 
26{xi + ai)e{x2 + a2)9[{xi - X2)9ixi + X2) , 

620 = - A6'(a;i + ai)9{x2 + a2)9{xi - X2)9{xi + X2) + ■ . . 

Proposition 6.2. For any Ki,K2 and generic 71,72 there is a unique (up to a 
constant factor) function $(a;) of the form (I6.6|l with the properties ()6.1()|l - (|6.13|l . 
It is described by the formulas H6.15|l - H6.18|l above. 

To prove the formula 16.15|l . one goes to the limit w in the formula (|7.7|l 
from the next section, picking up the first nonzero term (of order 6 in a>). 

Formulas H6.15|l - H6.18|l fix the dependence of $(x) on 4 parameters ai, 02, fci, ^2. 
Let us consider now the translation properties of 4> regarded as a function of these 
parameters. Recall that for generic ai,a2 the function $ was determined uniquely 
up to a factor by H6.6|l and (|6.12|) - H6.13|) . Thus, it follows immediately from Re- 
mark I^Jtha-t under the shifts H6.8|l - (|6.9|l <i> must remain the same, up to a factor 
independent on x. To find this factor, it is sufficient to look at the formula (|6.19|) 
for the leading coefficient &22- As a result, we conclude that the function $, given 
by formulas (|6.6|l . H6.15|l - (|6.18|l is invariant with respect to the shifts 16.8|l - (|6.9|l : 

(6.20) $(01 + 1, 02, ki, k2) = *(ai, 02 + 1, ki, k2) = $(ai, 02, fci, ^12) , 

(6.21) $(ai +T,a2,ki + 2m, k2) = $(ai, 02 + r, fci, fc2 + 27ri) = $(ai, 02, A:2) • 

Now let us find the interrelations between the parameters Oj , kj which will guar- 
antee two remaining vanishing conditions 1)6.10(1 - 16.1111 . Let Gi and G2 denote the 
derivatives 8182^ and 8182^ evaluated at xi = X2 =0: 

Thus defined Gi, G2 will be regarded as functions of the parameters ai, 02, ki, k2. 
Consider now the following two equations on these 4 parameters: 

(6.23) Gi=0, G2 = 0. 
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It is clear that conditions Ijti.lOII imply both of the equations (|6.23f) . Indeed, they 
guarantee that the function 

(6.24) /(i) =ai<f(0,i) 

is identically zero, in particular, /'(O) = Gi = and /'"(O) = G2 = 0. 

More interestingly, H6.23(l are 'almost' equivalent to H6.1()|l . To see this, note that 
the function H6.24|l has the following translation properties in t: 

f{t + l) = e^-}{t), 

From ^J^-^J^ we know that 

(9i±a2)$(0,0) = 0. 

This gives that /(O) = 9i$(0,0) = 0, while the first equation in (|6.23() gives that 
/'(O) = 0. Together with the translation properties above this implies that / is 
proportional to the following theta function: 

e''^'0{t + a2\TMt\T)f . 

This function has nonzero third derivative at i = as soon as 

(6.25) /c20(a2|r)+0'(a2|T) ^0. 

Thus, the conditions /(O) = /'(O) — /'"(O) — imply that / is identically zero 
provided that H6.25|l is true. The outcome is the following: for all ai,a2,fci,A:2 
satisfying the condition H6.25|l the equations H6.23|l imply (|6.10|l . 

Our next remark is that for any ai, a2, ki, k2 we have the relation did2^{0, 0) = 
df 82^(0,0). This is due to the identity 

45^92$ - 45ia|$ = (di + d2f{di - 82)^ ~ (di - d2f{di + 82)^ 

where the right-hand side vanishes at xi = 2:2 = due to (|6.12|) - (|6.13|) . Thus, 
repeating the same arguments as above, we conclude that the equations ()6.23|l 
imply also H6.11|l . as soon as 

(6.26) kie{ai\T) + e' {ai\T) ^ Q . 

As a result, we see that under assumptions (|6.25|I - H6.26|I both vanishing condi- 
tions H6.10|l - (|6.11|l are equivalent to one system (|6.23|l . To remove the restrictions 
(|6.25|l - (|6.26|l . let us consider equations (|6.23|) in more details. 

First notice that the functions Gi , G2 share the same translation properties 
(|6.20|I ~ H6.21|I in (ai, 02, /C2) with $ (since differentiating in x doesn't affect these 
properties). Notice also that Gi,G2 are polynomials in fci,fc2. Let us pass from 
fci, k2 to another variables pi,P2 as follows: 

(6.27) pi = fci -f C(ai), P2=k2 + Cia2), 

where we used C = C(^l''') to denote the logarithmic derivative of 0{z): 

e{z\T) 

This is slightly different from the Weierstrass C-function. 

Clearly, Gi , G2 are still polynomials in pi , p2 with the coefficients depending on 
oi, 02. The translation properties H6.20|l - 16.21|l imply that the coefficients in these 
polynomials will be elliptic functions of oi and a2- In fact, one can write down 
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Gi,G2 quite explicitly. The following result follows from our calculations for the 
difference case from the next section. 

Proposition 6.3. The system Gi = G2 = is equivalent to the following system: 

pM + K'2P2 + a') = p2{pi + K'lPi + cn , 

(6.28) pi{pi + m'2Pi + looi + (5cr + i5(a)2)P2 + cr + loaao 

= pM + locip? + loci'p? + (5cr + i5(a)')pi + cr' + iockh , 

where Ci, C2 stand for C(ai) ^'^'^ C(fl2) "ff/ji^e ^^e prime denotes taking the derivative 
with respect to the corresponding variable. (Notice that p{z) = —Q' + const, so 
Q" = — p' and so on.) 

From the discussion below will follow that for generic ai, a2 the system (|6.28ll has 
a finite number of solutions (pi , p2 ) ■ Thus , we can think of l().28|l as a finite covering 
of the product £ x £ oi two copies of an elliptic curve 5 = = C/(Z + rZ). In 
fact, the only (oi, 02) where the fiber is infinite are those with ai — ±02 (mod 1, t). 
This corresponds to the following 'vertical' components of (I6.28|l : 

(6.29) {ai = 02 , Pi = P2} , {ai = -02 , Pi = -P2} • 
Another 'trivial' component is, obviously, 

(6.30) {pi = P2 = 0} . 

If we delete these three components from (|6.28|) . the remaining part will be, in fact, 
a 13-fold covering of £ x £. Since we deleted the component H6.3()(l . the conditions 
and, hence, l|On|l - l|mT|l are satisfied on the remaining part. Thus, 
we arrive at the following theorem. 

Theorem 6.4. LetC be the finite covering of the product of two (punctured) elliptic 
curves which is obtained from (I6.28() by deleting the components (|6.29|l and (|6.30|l . 
Then C is the Hermite-Bloch variety for the operator H6.1|l and a double-Bloch 
solution ^{x) corresponding to a point (01,02,^1,^2) in C is given by the formulas 
(lO) . and (lOTjl . 

We still have to explain why C is a 13-fold covering. To this end let us consider a 
family of plane rational curves 1^9 : — > of degree 5, depending on a parameter 
a E £ = C/Z + tZ and defined as follows: ii a E £ and u — (uq : ui) E then 

ip{a, u) = {ipo ■■ (fii : (P2) , where 
ipo = uiUq , 

ipi = uful + 3C'(a)MiUo + C"(a)"o ' 

^2^ul + 10C{a)ulul + (5C"'(a) + 15{C {a) f)uiut, + {(""{a) + lOC (aK" {a))ul . 

Then the solutions (^1,^2) of H6.28|l correspond to the intersection points of two 
curves Ci = <^(ai, C2 — ^{0.2, ■ ) from our family. Namely, if (f{ai,u) — ^(02, v) 
then pi ~ ui/uq and p2 — vi/vq obviously satisfy H6.28|l and vice versa, provided 
P11P2 7^ 0. We should, however, exclude from consideration points with pi,p2 = 00. 
Namely, all the curves from our family pass through (0:0:1) = <p(a,oo). It is 
a standard exercise in basic algebraic geometry to show that mult (Ci 0^2) at 
this point equals 12. In doing this local analysis, one immediately observes that 
the condition C(ai) = Ci'^2) is necessary for Ci and C2 to coincide near the point 
(0:0: 1). This proves that the covering H6.28|l is finite apart from the components 
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(|5.29|l . After that Bezout's theorem teUs us that the number of common points, 
apart from (0:0: 1), equals 5 x 5 — 12 = 13. 

Remark 6.5. Operator H().l|) is symmetric under the Weyl group W of B2, whose 
8 elements act by permuting coordinates and/or changing their signs. This in- 
duces the action of W on Bloch solutions and, therefore, on the Hermitc-Bloch 
variety (|6.28|) . This action, in terms of (fli, a2:PijP2)j is generated by involutions 
(ai,a2,Pi,P2) {a2,ai,P2,Pi) and {ai,a2,Pi,P2) -> (-ai, 03, -pi,_P2)- 

6.1. Algebraic integrability. The operator L (|6.1|l is completely integrable. Ac- 
cording to Theorem 14. II it has a commuting operator of order four, 

= dfd^ + ... . 

Using Proposition 15. II we see that our ^ is a, common eigenfunction for L, Li: 

(6.31) Lip = El/;, Liilj = Eiip, L21P = E2^l) . 

Here E,Ei are some functions of the parameters oi, 02, /ci, which, in principle, 
can be calculated explicitly (though we didn't have enough energy to perform such 
a calculation). 

In 01] it was shown that apart from Li, there is another operator 

L2 = dl - hdldl + . . . 

which commutes with L, Li (see j44j for the explicit expression for Li, L2). The ex- 
istence of a fifth order quantum integral L2 means in this case that the Schrodinger 
operator L is algebraically integrable. To see this directly, let us consider one more 
operator L3, obtained from L2 by interchanging xi and X2- Then one easily checks 
that the common eigenspace (I6.31f) of L and Li is 8-dimensional, and for generic 
E and Ei it is spanned by the doublc-Bloch solutions tp, constructed previously. 
On the other hand, by proposition 15. II each ^|J will be an eigenfunction of L2 and 
L3 as well. So, the only thing to check is that the eigenvalues E2, E3 separate all 8 
solutions of the system H6.31|l . This is enough to check in the trigonometric limit 
T — !■ +ioo, which is not difficult. 

Remark 6.6. We do not give the precise relation between the Hermite-Bloch variety 
and the spectral surface. To find such a relation, a careful analysis of the structure 
of the divisor at infinity is needed. Let us remark that in '45' two algebraic relations 
between the 4 operators L, . . . , L3 were calculated explicitly. Thus, they determine 
a 2-dimensional affine algebraic variety in C*. However, it is not isomorphic to the 
spectral surface as an affine variety (though they are birationally equivalent). This 
can be seen already in the trigonometric limit t — > +ioo, by using the information 
about the spectral variety from 10 . Namely, the results of ^01 imply that these 
four operators L, . . . , L3 are not enough to generate the whole commutative ring 
(which is isomorphic to the coordinate ring of the spectral surface). 

6.2. Spectrum of L. Throughout this subsection we assume that the parameter 
T is pure imaginary, so the potential u{x) of the Schrodinger operator H6.ll) is real- 
valued for X S . Its singularities is the family of lines 

(6.32) Xi e Z , X2 (z Z , Xi + X2 (z Z , Xi — X2 G Z . 

These lines cut into triangles and the spectral problem for L splits into separate 
spectral problems for each triangle. 
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Let 'il'{x) = il}{x]ai,a2,Pi,P2) be a double-Bloch eigenfunction for L, which 
corresponds to a point (ai, a2,pi,P2) of the surface (|6.28|l in accordance with the 
formulas (|6.5|) . (|6.15|) - (|6.18|l and Ht).27|l . Given such a -0, let us symmetrize it in 
the following way: 

(6.33) -^{x) = ^ ip{wx) , 

where W denotes the Weyl group for the system 32- To get a square- integrable 
eigenfunction, according to Proposition 15.61 one takes k G 2TriP, that is k = 
(iTrm, iTTTi) with integer m, n having the same parity. For such k the substitution of 
(|6.27|l into 16.28|l leads to a system of (transcendental) equations on ai,a2, and by 
solving it one eventually finds the corresponding eigenfunctions. Note that because 
of the invariance of the system H6.28(l (and of -0) under the shifts H6.8(I ~ H6.9|I . it 
is enough to look for the solutions (0.1,02) with lying inside the fundamental 
parallelogram with the vertices ±(1 + t)/2. Also, taking into account the Weyl 
group action, we can restrict ourselves to the dominant weights, i.e. < m < n. 

Now let us consider the trigonometric limit t ~f +ioo, then one can show that 
for any m,n the corresponding system H6.28|l will have the unique solution 01,02 
inside the fundamental parallelogram (at least, for sufficiently big r). Moreover, 
if one fixes k = {ki,k2) and takes then the limit r +00, then will go to the 
Baker- Akhiezer function 'ip{k,x) considered in j42| (this can be seen directly from 
the formulas for tp). Now, for the Baker- Akhiezer function 'ip{k,x) it is known (see 
Theorem 6.7 of 42 ) that the formula (|6.33|l will produces all the Jack polynomials 
if 

(6.34) k = 27ri(A + p) with A e P+ and p - ^ ( nin. + l)a . 

For others k G the symmetrized ^ will be zero. 

In our situation this means that the function ^! = ^rn,n defined by will 
be non-zero as soon as 

(6.35) n — 4 > m > 2 , with n = m (mod2) . 

For other k E 2TTiP^ it will be zero in the trigonometric limit. But the result of 
cited in Section l5 . II claims that the family of the eigenfunctions of L is analytic in 
p = e'^''^ and specializes to the Jack polynomials at p — 0. Hence, if '^m.n = in 
the trigonometric limit, it must be zero for all r identically. 

We conclude that the eigenfunctions '^m.n of L are labeled by m, n satysfying 
()6.35|l . In particular, the ground state corresponds to {m,ri) = (2,6). The con- 
structed solutions ^'m_„ will have second order zeros along the lines (|6.32l) and will 
be invariant under orthogonal reflections with respect to these lines. According 
to the results of [311, the resulting family is complete in L^(r)^ (see Section |0] 
above). 

Remark 6.7. For certain (oi, 02, fci, fc2) the corresponding Bloch solution ■0 has a 
nontrivial symmetry in x, which must be a subgroup of the Weyl group W . Our 
formulas for ■0 do not work directly for some of these cases, because of the presence 
of an extra component H6.29|l . Of particular interest are those of the points which 
correspond to the solutions which are double-(anti)periodic (in each of the variables 
xi,X2)- These can be viewed as multidimensional analogues of the classical Lame 
polynomials |43| . 
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7. Difference B2 case 



In this section we wih generahze the results above to the following difference 
version of the operator (|6.1|) : 

(7.1) L = flo + a+Tf" + a_r-2'^ + fe+T|" + h^T-"^"^ , 

where stands for a shift in by e, and the coefficients a±, 6± are 



a± 

h+ = 



9{xi ± uj)6{xi + X2)0{xi — X2) 
d{x2 T i^)0{xi + X2T 2uj)9{xi - X2 ± 2cj) 



c± 
d+ = 



8{x2 ± uj)9{xi + X2)9{xi ~ X2) 
while ao has the form oq = c+ + c_ + + rf_ with c± , d± given by 

9{2uj)0{xi ± 5uj)9{xi +X2T 2uj)9{xi ~X2T 2a;) 
' 9{'iu;)9ixi ± Lu)9{xi + X2)9{xi ~ X2) ' 
9{2uj)9{x2 ± 5uj)9{xi +X2T 2uj)0{xi -X2± 2uj) 
9{Aiu)9{x2 ±uj)9{xi + X2)9{xi - X2) 

In all formulas 9{z) — 9{z\t) is the odd Jacobi theta function 15.4|l . 

This is a very special case of the so-called BCn generalization of the quantum 
Ruijsenaars model A7 ■ In trigonometric case it has been introduced by Koorn- 
winder 02]. Elliptic version was first suggested by van Diejen and later ex- 
tended in |5()j ■ where its complete integrability has been proven. This also can be 
viewed as an elliptic generalization of one of the Macdonald operators [3] for B2- 
In what follows we assume that the parameter u is generic. Note that the operator 
(|6.1|l can be restored (up to a certain gauge) in the limit uj ^ 0. 

It is worth mentioning that coefficients of L are not periodic, so instead of 
(double-) Bloch solutions one should look for the eigenfunctions in a certain 9- 
functional space. Another way of putting it is to observe that L can be reduced to 
elliptic form using proper gauge. For instance, consider 

L^S-^oLoS, 5 = 9{xi)9{x2)9ixi + X2)9{xi - X2) . 

Then L will have elliptic coefficients, so we can look for its Bloch solutions ip{x). 
Correspondingly, ^ — Sip will be an eigenfunction for L and it will have translation 
properties similar to those of S. Abusing the language, below we refer to $ as a 
Bloch solution for L. 

7.1. Bloch solutions. We are going to construct eigenfunctions of L similar to 
the differential case above. Our ansatz for $ remains unchanged: 



'i/3 




{3xi+ji\3t)9 







(3x2 -I- 72|3r) , 



(7.2) ^ = eiip{KiXi + K2X2) ^io^ 

0<i,i<2 

An analogue of the vanishing conditions Ht).10|l - (|6.13|l is dictated by the singularities 
of L and is the following: 

(7.3) $(w, 2:2) =$(-0), X2) for all X2 , 

(7.4) (^{xi,uj) =^{xi,-uj) for all xi , 

(7.5) '^{xi + ijj,X2 + uj) =^{xi ~ ll!,X2 — to) for xi + X2 ~ , 

(7.6) ^{xi + U!,X2 — o)) =^{xi — u),X2 + i^) for xi — X2 = ■ 
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We are going to show that for a certain two-dimensional variety in the space of pa- 
rameters 7i, 72, -ftTi, K2 there is only one (up to a factor) such (f>. As a consequence, 
$ will be an eigenfunction of L, due to a natural analogue of Proposition l5.ll 

We start from conditions H7.5() . (|7.6|l . Using the formula H6.4|l and repeating 
the arguments used in case a; = 0, we obtain a linear system for Cy and can see 
that (for generic 71,72) it defines Cij uniquely, up to a factor. However, solving 
this system leads to a very cumbersome formula. Instead, let us define $ by the 
following formula: 

(7.7) $ = exp(A:ia;i + k2X2)Y^ hj9{xi + ai + iuj)e{x2 + 02 + juj)e'^^'^'+^''^^ , 

where the summation is taken over the following set of indices: 

(z, j) = (0, 4), (4, 0), (0, -4), (-4, 0), (2, 2), (2, -2), (-2, -2), (-2, 2), (0, 0) , 
and the coefficients hij — hij [x) look as follows: 

bij = PtjO (^i +X2- ^ (^^1 - 2:2 - —Y'^ 

with 

/3o4 = /340 = /3o,-4 = /3-4,o = muj)f 

022 = /32,-2 = /3-2,-2 = P-2,2 = -0(2c^)0(4w) 

(3oo = m^)f . 

Proposition 7.1. For any Ki,K2 and generic 71,72 there exists unique (up to a 
factor) function $ of the form l|7.2|l satisfying conditions (|7.5|l - 1)7. It is given by 
the formula 1)7. 7|l . where kj — Kj -\- ni and Oj = jj — (1 -t- r)/2 (j — 1,2). 

To prove the proposition one first checks that this <i> has the needed translation 
properties in xi,X2, then an elementary check shows that conditions H7.5(l - ()7.6(l are 
satisfied. 

Let us turn now to conditions H7.3|l - (|7.4|l . First, let us remark that (|7.5|l . (|7.6|) 
at 2; = (0, 0) imply that 

(7.8) ^{u;,uj) - ^{-uj,-iLj) = ^{uj,-uj) - $(-cj,cj) = 0. 
Introduce now 

Gl = $(w, LU) - (f>(cj, -w) - <f>(-W, Oj) + $(-ti), -w) , 

G2 = *(t^, -3w) - $(-w, -3w) - $(w, 3w) $(-w, 3w) , 
and consider the system 

(7.9) Gi=0, G2 = 0. 

Obviously, the vanishing condition H7.3|) implies (|7.9|) . Conversely, from the system 
1)7. 9|) we deduce immediately that the function 

/(t) = a>(co,t)-$(-^,i) 

satisfies the conditions /(w) = f{—^) and /(Sw) = /(—3a;). Together with ()7.8)l 
this gives that /(oj) = f{^^) = 0. Since / is a theta function of order 3, it must 
have the form f — 9{t — ijj)9{t + uj)g{t) with g{3uj) — g{~3uj). Now, since g is a 
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theta function of the first order with known characteristics (expressed in terms of 
a2, ^2), the condition g{5uj) — g{—3Lu) imphes g = as soon as 

(7.10) ^'^''^PrPn- 

The latter condition is a difference version of (|6.25|l . Now let us collect some 
corollaries of (|7.5() - (|7.6() : 

-3w) = $(3cj, ~Lu) , ~Slu) = $(-3cj, -LU) , 

3w) = $(3tJ, tj) , 3cj) = $(-3w, w) . 

Thus, G2 can be rewritten as 

G2 = «'(3w, -cj) - $(-3tJ, -cj) - $(3w, w) + $(-3w, . 

Hence, we can repeat the same arguments with respect to X2-variable and conclude 
that the system (|7.9|l implies (I7.4|l as soon as 



(7.11) e^-^^ ^ 

0[ai + Suj\t) 

Summing up, we see that the system (17.911 implies both of the conditions (|7.3|l - (|7.4(l 
under assumptions H7.10|I - H7.1HI . 

To get rid of restrictions H7.10|l - H7.11|l let us calculate Gi,G2- First, introduce 
the notation ^1,^2 for 

A direct substitution gives that 

$(w,w) ^PQ^-4:O{'iLo)9{~2uj)0{ai + uj)9{a2 ~ 3uj)^i^2^+ 

p-4,ae{4uj)e{2uj)e{ai - 3uj)9{a2+iu)^i^^2+ 

P-2,20{2uj)e{2u)e{ai - u)e{a2 + 3cj)er^C|+ 

(32,~20{2uj)e{^2w)e{ai + iuj)e{a2 - cj)£.U2^ ■ 

In a similar way we calculate ^{—u!,u!) and Gi — 2^{u!,u!) — 2$(— w). As a 
result, the equation Gi = takes the following form (up to a nonessential factor): 

(7.12) {e{ai + 3lj)£_1 - e{ai - 3Lj)(^^){0{a2 + 1^)^ - 0{a2 - lu)C2^) 

= {9{ai + uj)Ci - e{ai - uj)^^^){e{a2 + 3uj)^l - 6(02 - 3a')C2~^) . 

In a similar way one can compute G2 . It turns out that it is a linear combination 
of 8 terms of the form ^i^^2^ and 8 terms of the form 

Moreover, the combination of the last 8 terms is proportional to Gi, so we can get 
rid of them by subtracting Gi, this does not affect the system ()7.9|l . After these 
transformations equation G2 = takes the following nice form: 

(7.13) {e{ai + 5lj)£.1 - 0(ai - 5Lj)(^^){0{a2 + 3lu)£_^ ~ ^(02 - 3lu)(^^) 

Summarizing, we see that the system H7.9|l is equivalent to the equations H7.12|l - 
(I7.13|l . These equations are obviously invariant under the transformations (|6.8|l - 
(16. 9|) , in this way they define a covering over the product 5 x £ of two elliptic curves 
£ = C/Z + tZ. It has two 'vertical' components 

(7.14) {ai = 02 , a = 6} , {ai = -"2 , 6 = (6)"'} . 
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Another 'trivial' component is given by 

(7.15) 9{ai + 3uj)^l - e{ai - 3w)Cr^ = 6l(a2 + 3uj)(l - e{a2 - 3w)C2~^ = . 

After deleting these three components, one gets a finite (in fact, 17-fold) covering 
of £ X £, let us denote it by C We can conclude now that for any point in C 
the corresponding function <I> will satisfy the vanishing conditions (|7.3|l - l|7.6|l and, 
hence, it will be an eigenfunction of the difference operator L. 

Theorem 7.2. Formulas H7.7|) and equations (|7.12|l - (|7.13|) describe the (double-) 
Block eigenfunctions of the difference operator l|7.1|l . The Bloch-Hermite variety 
C, which is obtained by deleting components H7.14|l - 117. 15() from the variety (|7.12ll - 
(|7.13|l . is a 17 -fold covering of the product £ x £ of two elliptic curves. 

As a corollary, considering the limit lu ^ we can calculate explicitly the 
Hermite-Bloch variety for the operator H().l|) . Namely, one picks up the terms 
of order 4 and 6 in w in equation (|7.12l) . 

Corollary 7.3. In the limit uj the system H7.12(l - (|7.13|l goes to H6.28(l . 

The only thing we still have to explain is why the degree of the covering is 17. 
To this end, let us define a family of plane rational curves </? : — > of degree 5, 
depending on parameter a G £. Namely, for a ^ £ and u = {uq : ui) G P^ put 

ip{a, u) = {ipQ : ipi : (^2) , where 

ulul0{a + uj)0{a - uj/2) ulul0{a - u!)0{a + uj/2) 



= 



9{a + uj/2) 0{a-w/2) 

uouj0{a + 3uj)0^{a - uj/2) uffUie{a - 3uj)e^{a + uj/2) 

e^{a + uj/2) 0^{a-uj/2) 

u^0{a + 5uj)0^{a - uj/2) u^6l(a - 5uj)0^ia + uj/2) 



95{a + uj/2) 0^{a-uj/2) 

Then the solutions (^1,^2) of (|6.28|) correspond to the intersection points of two 
curves Ci = ip{ai, ■ ), C2 — '^{0,2, ■ ) from our family. Namely, if ip{ai,u) = ip{a2, v) 
then ^1, ^2 with 

ui0'^{ai - Lu/2) 2 vi0'^{a2 - Lu/2) 



iiir - , and (e 



uo0^iai+Lu/2) vo0Ha2+Lo/2) 

clearly satisfy H7.12|l - (|7.13|l and vice versa, provided H7.10ll - H7.11|l . We should, 
however, exclude from consideration points with ^1,^2 — 0, 00 since — e"*^'. 
Namely, all the curves from our family pass through (0:0:1) = (p(a, 0) = ^{a, 00). 
A little difference with uj — case is that now we have mult (Ci fl C2) = 8 at this 
point. By Bezout's theorem, the number of intersection points of Ci,C'2 , apart 
from (0:0: 1), equals 5 x 5 - 8 = 17. 

Remark 7.4. Note that for given 1^^, ^| the corresponding quasimomenta ki, k2 seem 
to be non-unique, with the ambiguity of adding some multiple of ii: /uj. However, 
this would result in multiplying $ by e'^^^^/'^ ^^n^x^/uj ^j-g quasi-constant on 

the lattice 2ujl?. Thus, this leads to the same eigenfunction, so the constructed 
Bloch solutions are in one-to-one correspondence with the points of the surface C. 

Remark 7.5. It is clear that the Weyl group action on C is generated by two invo- 
lutions 

(ai,6) ^ (-ai,Cr^) and (ai,^i) ^ (02,6) • 
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7.2. Structure of the solution space. As we mentioned above, the difference 
operator (|7.1|) is completely integrable, i.e. there exists another difference operator 
Li which commutes with L. It is given by the following expression |49[ I5U) : 

(7.16) Li = c++T^T!^ + c+^T^T^"^ + c^+T^'^T^ + c^^T^'^T^'^ , 
where the coefficients Ccj.c2 look as follows (we treat ± as ±1): 

_ Ojxi - ^61)^(3:2 - Lue2)0{xi + X2- 2uj{ei + e2))d{xi ~ X2 ~ 2tj(ei - 62)) 

^'^''^ " e{xi)e{x2)e{xi + X2)e{xi - X2) ■ 

Now let us consider the system of two partial difference equations: 

(7.17) Lf^Ef, LJ^Eif, 
defined on the lattice 2u!£ where 

C = {{m, n) I m ± ri G Z} . 

More precisely, we fix generic x^ g as a base point and regard a function / in 
(|7.17|l as being defined on x^ + 2bjC C . The base point must be outside the 
singular locus of L, Li, i.e. such that L, Li are nonsingular on x^ + 2ujC. 

The Bloch solutions constructed above, are common eigenfunctions of L and 
Li. (The proof for Li is the same: the only thing to check is an analogue of 
Proposition 15. II ) Since L and Li are ly-symmetric, each of the 8 functions ^{wx) 
{w S W) will solve the system (I7.17|l . We know that for generic point of the spectral 
surface X (thus, for generic E,Ei) all 8 functions ^{wx) are linearly independent 
(as functions on C^), because they have different translation properties with respect 
to the shifts Ht).8|) - (|6.9|) . Hence, their restriction to a;° + 2ujC also gives 8 linearly 
independent solutions of (|7.17() fat least, for generic base point a;°). On the other 
hand, it is not difficult to see that any solution / is uniquely determined by its 
values at eight points x'^ + v with the following v. 

V - (0, 0) , (±^, ^) , (lo, ^lu) , (±2w, 0) , (0, 2lo) , {to, 3lo) . 

This implies that the dimension of the solution space of l|7.17ll is at most 8. Thus, 
we conclude that any solution of (|7.17(l (for generic E, Ei) is a linear combination 
of 8 Bloch solutions {^{wx)}^^w 

Proposition 7.6. The space of solutions of the system H7.17|l has dimension 8 and 
for generic E,Ei is generated by the Bloch solutions ^(wx) (w € W). 

Remark 7.7. Above we associated a double-Bloch eigenfunction $ to a solution 
(fli, a2, ^1, ^2) of the equations H7.12|l - H7.13|) . Note that these equations are invari- 
ant under —^j, but this does not lead to another eigenfunction, since they will 
differ by a quasiconstant factor. Situation is different for the system (17.17(1 . since 
it is defined on a different lattice. It is easy to see that (^1,^2) and (— Ci, —£.2) still 
lead to the same solution modulo quasiconstants, the same is true for (^1, —^2) and 
(—^1)^2)- The resulting two functions have the same eigenvalue E in H7.17|l . but 
opposite values of Ei. Thus, the Bloch variety for the system (|7.17|) is a double 
covering of the surface C introduced above. 

We will not go into discussing the spectral properties of the difference operator 
L. See papers [521 ESI devoted to this rather delicate matter. Let us just remark 
on some special solutions analogous to the 'discrete spectrum' considered in section 
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15.21 Namely, let us consider the following anti-invariant solution of the system 
(TTTTIi : 

^skew{x) = ^ {det w)^{wx) . 

The vanishing conditions (|7.3|I - H7.6|I imply that ^skew vanishes along lines xi ± 
X2 — ±2uj and Xj — ±a;. It also vanishes if xi = ±X2 and xj = due to anti- 
invariance. Let us require now for all 8 functions ^{wx) to have the same Floquet- 
Bloch multipliers with respect to the shifts by ei and 62, which is equivalent to the 
conditions 

exp(fci) = exp(fc2) — ±1 • 

Then ^skew will vanish also along the shifted lines 

Xi ± X2 = m , m ± 2u! , Xj = n ,n ± to [m, n € Z) . 

In the limit w ^ these solutions go to those 4* constructed in section more 
precisely, 

oJ-^-^skew — > M{9'{Q)fe{ai)9{a2)e{xi)e{x2)9{xi + X2)e{xi - xa)^' as oj ^ . 



8. HiETARINTA OPERATOR AND ITS DISCRETIZATION 

8.1. Continuous case. We consider now the Schrodinger operator (|1.4|l but first 
let us rescale the coordinates Xi — > aiXi, so instead of l|1.4|) we will consider 

(8.1) L = -ajdf - aldl - aldl 

+ 2{a{ + a\)p{x\ - X2) + 2(a2 + a3)p(a;2 - xz) + 2(a3 + a^^p[xz - x\) , 

where, as before, p{z) = p(z|t) is the Weierstrass p-function and -|-a| +03 = 0. 
We are going to calculate the double-Bloch eigenfunctions of L. More specifically, 
we are looking for the solutions il) of the equation L'0 — Etp with the following 
properties: 

(i) ij; is of the form 

^(x) 

(8.2) Tpix) = — — — exp(A:iXi + fc2a;2 + fc3X3) , 

6'(a;i2)6'(a;23)6'(a;3i) 

where 6 — 9^lJ'^ and $ is holomorphic in C"^ and depends on the 

differences x^ only, in other words, {d\ + ^2 + ^3)$ = ; 

(ii) i\} has the following translation properties: 

(8.3) ^{x^'ej)^e^^i){x), ^{x + Tej) ^ e^'^ i)(x) (j = l,2,3), 

where (ei, 62, 63) is the standard basis in C^. 

It is not difficult to conclude that for fixed kj , puj the conditions above determine 
a three-dimensional functional space, and the corresponding (f>(x) in H8.2|l must be 
of the form 

2 

(8.4) * = ^ Q^(a;i2 + 612 + /t/3)0(x23 -I- 623 + ?r/3)0(a;3i + &31 + ^r/3) , 
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where co,ci,C2 are arbitrary constants and the parameters &12, ^23,^31 are related 
to /ij above in the following way: 

gM2 _ pfe2r+27rjbi2-27r'i623 
gM3 _ gfe3r+27rjti23-27r'i&3i 

This shows that the three-dimensional space (18 .41) depends, essentially, on the pair- 
wise differences of the parameters bim only. Thus, without loss of generality we 
may assume that 

(8.6) 612 + &23 + &31 =0. 

In formulas below we will also use 621, ^32, ^13 under the convention that 6^ = ^^ji- 
Now, in accordance with Proposition l5.ll we impose certain vanishing conditions 
on ijj which are motivated by the structure of the singularities of the operator 
(|8.1|) . Namely, for any i — 1, 2, 3 consider the function f{t) — ■ip{x + taf_^ei^i — 
taf+iCi+i) for X such that Xi^i — Xi+i (we treat indices modulo 3, so xq = 2:3). 
Our assumptions about -0 imply that for such x the function / will have a pole 
at f = 0, so its Laurent expansion will look as / = a-it~^ -I- ao -I- aii -I- .... The 
coefficients in this expansion depend on x. Let us require that ao = for all x such 
that Xi-i = Xi+i. Using (|8.2I) one rewrites this condition as follows: 

(8.7) af^^ki^i - af^-^ki+i + = , 

p _ af_i9,-i$-af+i9,+i$ ^ 6>'(x,_i,,) 3 0'{x.,+i,^) 

with l|8.7|l to be valid for all x such that Xi-i — x^+i. 
The following lemma follows from Proposition 15. II 

Lemma 8.1. If ^ has the form H8.2() . (|8.4|l and satisfies the vanishing conditions 
(|8.7I) . then the same will be true for its image t/j = Ltp under the action of the 
operator (|8.1|1 . 

As we will see below, for a certain three-dimensional subvariety in the space of 
the parameters kj , 5/^ the vanishing conditions cut a one-dimensional subspace in 
the space (|8.4|l . Thus, the lemma ensures that the corresponding ■ip{x) will be an 
eigenfunction of L. 

We may regard the restriction of the expression Fi on the plane Xi-i — Xi+i as 
a function of z = Xi-i^i = Xi+i^i. It is easy to check then that Fi will be an elliptic 
function of z with periods l,r. So, first of all we have to choose the parameters 
bim,Cj in such a way that Fi{z) would be non-singular. Let us assume that the 
parameters 01,02,03 are generic enough, i.e. that of ^ aj. Then for Fi to be 
non-singular at z = we need $(0) = 0. This gives the following condition: 

(8.9) CQ + ci + C2 = 0, ci^cie{bi2 + ^-^)0{h2z + %)e{hi + %). 



Now we note that 



i+i 



is a one- dimensional 0- function of order 2, hence it has two zeros (modulo 1,t). 
First zero is z = (due to condition H8.9|l ). An easy check shows that the second 
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zero is z = + So, up to a constant factor, 

$(z) = 0(z)0(z-6,^,_i-6,,,+i). 

To get rid of a possible pole at z = + in (|8.8I) we must require that 

af_idi-i^ — af^j^di+i^ = for Xi-i,i = Xi+i^i = + This leads to the 

following relation: 

2 

(8.10) J2 M<^lC{b23 + alCibsi + y ) + alC(&i2 + |)) = • 
;=o 

Here and below (^{z) ■j^log6{z). Notice that the relation H8.10|l is symmetric 
with respect to indices 1,2,3 (so we have just one condition instead of possible 
three!). 

We use the relations (|8.9|) . (|8.1U|) to express (up to a common factor) the param- 
eters ci in terms of the parameters 5^ . These relations imply that each of Fi (|8.8I) is 
nonsingular in z Xi-i^i ~ Xi+i.i^ therefore they are some constants depending on 
ci,bij. Thus, H8.7(l leads to the expressions for the differences a'^_iki-i — af_^_iki+i in 
terms of c/ and bim- One can check that the resulting system is always compatible 
(i.e. that Fi + F2 + = 0). This follows, for instance, from the compatibility of 
the system (I8.18|) below by going to the limit — > 0. As a corollary, the formu- 
las (|8.2() . H8.4|I . together with (|8.9|l . H8.10|l and (|8.7() deliver the expression for the 
double-Bloch eigenfunctions of the operator L. 

Finally, let us discuss the structure of the Hermite-Bloch variety of the op- 
erator H8.1|l . The double-Bloch solutions are parametrized by foi2,&23,^3i with 
bi2 + 623 -I- 631 — 0, and the corresponding ki, ^2, fca are determined from H8.7|l . De- 
note by b and k the three-component vectors b = (612, 623, 631) and k — (ki, fc2, k^). 
Then two different points in the parameter space (5, k) lead to the same solution 
iff the corresponding Floquet multipliers in H8.3|l are the same. Taking into ac- 
count relations H8.5(l . we conclude that the following transformations do not lead 
to another Bloch solution: 

b^b + si, k-^k (£1 (2/3,-1/3,-1/3)) 

(8.11) 5^6 + 62, k^k (£2 = (-1/3,2/3,-1/3)) 
5^5-|-T£i, fc-f 27ri(l,-l,0) 

5^6 + T£2, A: -> fc + 27ri(0, 1, -1) . 

Thus, b is effectively represented by a point of a factor C'^/ZI + tC where €? = 
{zi + Z2 + Z3 = 0} C and the lattice C is generated by £i,£2. This factor is 
isomorphic to the product of two elliptic curves with parameter r. Above each 
point b we have a complex line of double-Bloch solutions, because equations H8.7|l 
determine k up to adding any multiple of (a|j~^, a^^, a^^)- 

8.2. Discrete case. We keep the notation Xij for Xi — Xj. The discrete version of 
the operator l|8.1() looks as follows: 

3 



■ ^ ^ e{cjaf)e{x,^,,,)9{x,,,+,) 



where al + + — and stands for a shift by e in a;^. Its rational version 
9{z) = z was communicated to us by M.Feigin who found it to be dual (in bispectral 
sense) to the trigonometric version p — sin^^ of the Hietarinta operator. 
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The difference operator (|8.12|) relates to (|8.1|) in the foUowing way: 



13 = ^ + a2 ^ + a3 ^ + uj{di + ^2 + d^) + ^(const ~ L) + o{uj'^) as a; ^ , 

where L is gauge-equivalent to L, 

L = 5 o L o S^^ , S = 6{xi — X2)9{x2 — xj,)9{xz — xi) . 

Unlike L, the operator D is not periodic. As a result, instead of the double 
Bloch eigenfunctions, we will look for eigenfunctions with the translation properties 
similar to those of 5. Apart from that, our ansatz for the eigenfunctions of the 
operator D remains the same: 

(8.13) = exp(fcixi + ^2X2 + fcaza)^ , 

2 

(8.14) * = 51 ^(^(2:12 + + /t/3)0(x23 + ^23 + ^t/3)0(x3i + ^31 + Ir/i) , 

(8.15) 612+623 + 631=0. 

The vanishing conditions now look as follows: for each i = 1, 2, 3 

(8.16) := e{x,,,^^+ujati)T^i^\^)-e{x,,,+,+u:aj+,)T^^\^) = 

identically for all x with Xij^i — Xi^i. 

We have then a straightforward analog of Lemma [8. II so the same approach as 
above will give us the eigenfunctions for D. 

Let us first formulate the result. Namely, we consider the following two condi- 
tions on the function ^{xi,X2:Xz) given by (|8.14(l : 

(8.17) $(tja?,0,-cja^) = 0, <^{~ujal,{),ujal) = . 

This gives us two linear equations on c; and we use them to express c/ (up to a 
factor) through 6^. 

Secondly, we impose the following three relations: 



2 



e 



>alk,-u.alk, _ g(^ai)$(0,0,^ai) 
e{uja{) $(waf,0,0) 



;aifc3— a^fc. _ ^(wai) $(0, Wa^, 0) 



e 



We use these formulas to express ki,k2, through <&. The solution is not unique, 
and in fact we have a one-parameter family of fcj. Altogether, formulas H8.17|) - l|8.18|) 
fix the dependence of c/ and ki and hence of tp on three parameters bij (related by 
(|8.15|l '). The resulting family of functions ip{x) depends on three parameters: two 
of bij and one more due to the freedom in resolving (|8.18l) . see more comments 
below. 

Theorem 8.2. The formulas and the relations l|HT7|l - (|^TTS|l give a 

three-parameter family of eigenfunctions for the difference operator D. 

To prove the theorem, let us first notice that each of Fi in H8.16I) , being regarded 
as a function of z = Xi^i-i = Xi^i^i, is a one-dimensional theta-function of order 3, 
so if it doesn't vanish, it must have three zeros (modulo l,r). Moreover, a simple 

35 



count shows that the sum of these zeros will be equal to + On the other 

hand, a direct substitution into H8.16|l shows that the relations H8.17|l imply that F2 
vanishes for z = —toal and z = —loo^. Further, the first relation in H8.f 8|l simply 
encodes the fact that F2 vanishes at z = 0. Since the sum of these three zeros is 
uja2 which, generically, is not 621 + ^23, we conclude that F2{z) is zero identically. 

At first glance it seems that we need to add four more conditions to ensure all 
the vanishing properties (|8.16() . Namely, one needs also 

(8.19) $(0, ival, -ujaj) = $(0, -ujaj, ujaj) 

= $(-^03, ual, 0) = ^{toal, -cjflg, 0) . 
However, since $ depends on the pairwise differences of Xi only, we will have that 
$(0, 0702, —ujai) — wflj, —Lual + uoa\, ~oja\ — LLia\) = (^{—ULia\, 0, uoa^) — . 

In the same way other relations in (|8.19() follow from (|8.17|) . 

This demonstrates that the three-parameter family constructed in the theorem 
satisfies the vanishing conditions H8.16|l . thus proving the theorem. 

Finally, let us comment on the structure of the Hermitc-Bloch variety. Similarly 
to the case uj — Q above, Bloch solutions are parametrized by 6 = (612, 623, 631) with 
612 + ^^23 + ^31 = 0. This determines the corresponding c/ by H8.17|l . After that 
k — {ki, k2, ks) are determined from (|8.18l) . At this point we have certain freedom: 
if k = {ki, fe, fcs) is a solution of H8.18|l . then any 

, ( , t , 2 2 2\ StTZ , 2 2 2\ 

k = k -\ [a^ , a2 ,03 jH (nia-^ ,71202 ,71303 j 

UJ CO 

with any t (z C and integer ni, 712, 713 will be a solution, too. However, the last term 
is not essential since it results in multiplying $ by a quasiconstant. For the same 
reason, the factor t in the second term is essential modulo 27rj only. Besides, we 
still have the translation invariance of $ with respect to the transformations H8.11|l . 
Thus, the Hermite-Bloch variety is fibered over the product of two elliptic curves 
with the fibers isomorphic to C/27ri. 
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